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ABSTRACT 

Two  computer  methods  for  the  design  of  radial  out-flow 
turning  vanes  are  described.   The  key  requirement  was  for 
vanes  with  simple  shapes  that  could  be  machined  easily  and 
inexpensively  in-house.   Two  design  approaches,  one  using 
circular  arcs  for  both  suction  and  pressure  sides  and  one  for 
wedge-arc  shaped  sections  were  developed. 

Requirements  for  the  general  design  procedure  arose  in  two 
exploratory  projects,  a  Rotary  Detonative  Wave  Engine  <RDWE) 
and  a  high  speed  Centrifugal  Diff user  Test  Device  (CDTD) .  The 
RDWE  was  analyzed  first  and  found  to  be  potentially 
inefficient  in  its  original  configuration.  Preliminary  design 
for  a  high  speed  CDTD  flow  generator  was  completed, 
incorporating  wedge-arc  radial  outflow  vanes.   An  attempt  to 
analyze  the  flow  through  the  geometry  using  the  Eidelman- 
Godnov  Euler  (EGE)  was  initiated. 
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NOMENCLATURE 

A»  throat  area 

Ah  throat  area  for  High  Speed  CDTD 

Ai  throat  area  for  Low  Speed  CDTD 

a  throat  dimension 

b  width  test  section 

D  outer  diameter 

e  energy  of  the  unit  volume 

L  length  of  blades 

m  mass  flow  rate 

P  power 

Ps  specific  power 

p  pressure 

R  radius  of  curvature 

ri  inner  radius 

r2  outer  radius 

s  blade  spacing 

T  torque 

U  velocity  component  corresponding  to  x  coordinate 

Un  velocity  component  normal  to  edge 

U2  wheel  rim  velocity 

V  velocity  component  corresponding  to  y  coordinate 

Vfc  velocity  component  tangent  to  edge 
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V  9i  tangential  component  of  inlet  velocity 

V ©2  tangential  component  of  outlet  velocity 

W2   relative  exit  velocity 

Wx    work 

w'   rotational  velocity 

Z    number  of  blades 

a    efflux  angle 

82   blade  exit  angle  with  respect  to  radial  direction 

B2'  negative  blade  exit  angle  with  respect  to  radial 

direction 
T    specific  heat  ratio 
£    internal  energy 
"c   power  conversion  efficiency 
£   density 

Notations  used  in  Appendices  B  and  C  are  independently 
defined  therein. 
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I.   INTRODUCTION 

The  motivation  for  the  work  reported  here  came  from  two 
exploratory  projects,  a  Rotary  Detonative  Wave  Engine, 
(RDWE),  described  in  Ref .  1  and  a  high  speed  Centrifugal 
Diffuser  Test  Device  <CDTD),  described  in  Ref.  2.   Both 
projects  required  the  geometrical  design  of  an  axisymmetric 
generator  for  nearly-tangential  transonic  outflow. 

To  improve  gas  turbine  engine  performance  an  increase 
in  cycle  pressure  ratio  is  needed.   One  alternative  to 
mechanical  compression  is  to  supply  the  compression 
thermally  by  detonative  combustion.   If  the  pressure  rise 
produced  in  this  manner  is  sufficient  for  the  engine,  no 
compressor  and  associated  drive  turbine  are  required  and  a 
simpler  more  efficient  engine  might  result. 

As  described  by  Honks  [Ref.  13  ,  a  proposed  geometry 
for  such  an  engine  might  include  a  rotor  with  outflow 
passages  which  are  turned  at  the  exit  to  produce  nearly 
tangential  flow.   Figure  1  illustrates  the  proposed  RDWE. 
Honks  attempted  to  measure  the  torque  on  a  static  simulated 
rotor  as  a  detonation  wave  propagated  from  the  center 
outward.  A  reaction  torque  was  produced  from  the  nearly- 
tangential  expulsion  of  gases  from  the  rim. 
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The  question  then  was,  what  would  be  the  best  design 
for  the  rotor  passage  to  maximize  the  torque  impulse?  In  an 
analysis  included  in  Appendix  A  it  was  shown  that  tangential 
flow  at  the  rim  of  the  RDWE  would  yield  the  highest  specific 
power  and  highest  power  conversion  efficiency.  A  design  was 
required  therefore  of  nozzles  which  would  direct  the  outflow 
almost  tangent ially .  The  nozzles  would  operate  in  an 
unsteady  mode. 

The  second  project  which  motivated  the  present  study 
involved  the  development  of  a  test  rig  for  radial 
diffusera.   Such  diffusers  are  used  in  centrifugal 
compressors  to  convert  the  kinetic  energy  of  the  flow  from 
the  rotor  into  static  pressure.   Erwin  CRef .  2]  proposed  and 
designed  a  so-called  low  speed  Centrifugal  Diffuser  Test 
Device,  (CDTD) .   The  CDTD  simulated  the  flow  of  air 
delivered  to  the  diffuser  from  the  compressor  rotor  using 
static  swirling  vanes.   The  facility  would  allow  the 
evaluation  of  proposed  new  diffuser  geometries  and  permit 
validation  of  computer  analyses  for  diffusers  operating  in 
the  fully  subsonic  flow  regime.   Figure  2  shows  a  schematic 
of  the  low  speed  CDTD.   Vidos  CRef.  3]  presented  results 
from  the  flow  generator  used  in  the  CDTD  and  performed  an 
analysis  of  the  internal  flow.   As  mentioned  by  Vidos,  the 
design  of  centrifugal  diffusers  is  presently  based  on 
experimental  results  for  two  dimensional  and  conical 
diffusers.   Numerical  methods  are  currently  under 
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development  for  viscous,  three  dimensional,  transonic  flow 
(a  regime  containing  both  supersonic  and  subsonic  regions) , 
with  adverse  static  pressure  gradients.   Such  methods  are 
required  to  design  optimum  radial  diffusers  for  transonic 
conditions.   Equally,  a  high  speed  version  of  the  CDTD  was 
required  in  order  to  examine  such  diffuser  flows  under 
controlled  conditions  and  to  verify  the  emerging  codes. 
Again,  nearly  tangential  swirling  outflow  was  required,  this 
time  in  steady  flow. 

The  main  interest  here  was  therefore  to  design  vane 
shapes  and  the  apparatus  for  a  high  speed  CDTD.   The  design 
would  closely  parallel  the  general  arrangement  of  the  low 
speed  CDTD.   The  high  speed  CDTD,  had  the  unusual 
requirement  of  radial-outflow  with  convergency  at  the  rim 
exit.   The  limitation  on  the  design  was  primarily  that  the 
hardware  could  be  machined  easily  and  inexpensively  in- 
house .   Computer  programs  for  the  design  of  supersonic 
nozzles  for  radial-inflow  turbomachines  were  found  to  exist, 
but  most  were  highly  involved  and  would  need  extensive 
machining  processes  to  generate  the  blade  shapes. 
Particular  designs  are  discussed  herein.   A  literature 
search  indicated  that  computer  programs  for  radial-outflow 
devices  using  simple  vanes  for  two-dimensional  isentropic 
flow  of  a  perfect  gas  were  not  available. 

In  the  present  report,  the  design  of  the  blading  for  a 
high  speed  CDTD  model  and  the  general  arrangement  of  the 
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apparatus  are  described.   General  methods  and  computer 
programs  for  the  design  of  double  circular  arc  and  wedge  arc 
blading  for  radial-outflow  vanes  are  described  in  detail  in 
appendices.   The  report  also  discusses  a  method  of 
predicting  the  flowfield  through  a  selected  passage  using  an 
Euler  code  based  on  the  Godunov  method  of  solution. 
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II.  APPARATUS  DESIGN 

Factors  considered  in  the  design  of  the  high  speed  CDTD 
fell  into  two  categories;  namely,  limitations  due  to  the 
available  air  supply  system  and  piping  and  the  requirement 
that  the  test  section  have  dimensions  appropriately  scaled 
for  the  diff users  required  to  be  tested. 

A.   DESIGN  CONSIDERATIONS 

The  high  speed  CDTD  was  designed  to  operate  as  an 
intermittent  blow-down  wind  tunnel.   In  order  to  use  the  air 
supply  system  available  in  building  230  at  the  Naval 
Postgraduate  School,  the  test  section  flow  area  was  dictated 
by  the  available  mass  of  stored  air  and  the  minimum  desired 
runtime  of  two  minutes.   Using  proportions  similar  to  those 
used  in  the  low  speed  CDTD  (see  Figure  2),  the  radius  of  the 
swirl  generator  became  four  inches  while  the  radius  at  the 
test  section  was  five  inches.   The  method  used  for 
calculating  these  dimensions  is  given  in  the  following 
section.   The  difference  between  the  swirl  generator  outer 
radius  and  the  test  section  radius  allowed  for  mixing  of  the 
wakes  from  the  blading  (or  nozzles)  used  in  the  swirl 
generator.   The  swirl  generator  was  limited  to  a  length  of 
two  inches  to  meet  the  scaling  requirements  set  by  the  test 
section.   Because  of  the  larger  pressure  differential 
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required  across  the  swirl  generator  in  the  high  speed 
version,  turbine  like  turning  vanes  were  selected  rather 
than  the  sheet  metal  nozzle  construction  used  in  the  low 
speed  device. 

B.   APPARATUS  DESCRIPTION 

Figure  3  illustrates  the  overall  arrangement  of  the 
proposed  CDTD  design.   The  swirl  vanes  for  the  high  speed 
CDTD  are  sandwiched  between  two  circular  plexiglas  disks. 
Plexiglas  disks  were  used  to  allow  visualization  of  the  flow 
through  the  vanes.   The  disks  are  mounted  on  an  aluminum 
center  pipe  shown  in  Figure  3.   The  pipe  provides  double 
axial  in-flow.   Following  the  work  performed  by  Feieresien 
[Ref .  4]  a  transition  section  to  turn  the  air  from  the  axial 
to  the  radial  direction  was  added  in  the  form  of  a  plate  on 
a  tie-rod  on  the  pipe  centerline.   The  design  of  the 
transition  section  is  not  complex  if  the  area  of  the  radial 
outlet  is  smaller  than  that  of  the  total  axial  inlet.   An 
acceleration  through  the  transition  section  will  then  occur, 
leading  to  generally  favorable  pressure  gradients.   These 
consideration  required  a  double-entry  supply  arrangement  for 
the  CDTD  model. 

The  piping  is  attached  to  the  air  supply  system.   Figure 
4  shows  the  configuration  of  the  blades  on  the  plexiglas 
disks  and  in  relationship  to  the  supply  pipe. 
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C.   DESIGN  CALCULATIONS 

Vavra  CRef .  5]  showed  that  in  the  design  of  a  turbine 
cascade  the  efflux  angle  a  can  be  calculated  from  the 
geometry  of  the  cascade  and  that  it  is  independent  of  the 
incidence  angle  i  of  the  inlet  velocity,  see  Figure  5. 
Vavra  also  noted  that  it  is  common  practice  in  steam  turbine 
design  to  use  the  relation 

a  =  arc  cos  <  a  /  s)  <1) 

where  "a"  is  the  throat  opening  at  the  discharge  and  s  is 
the  blade  spacing.  Equation  (1)  applies  when  the  trailing 
edges  of  the  blades  are  of  negligible  thickness.   In 
practice  the  "wedge"  angle  at  the  trailing  edge  and  trailing 
edge  bluntness  modifies  this  dependence. 

With  the  limitation  of  a  two  minute  runtime  and  the 
available  mass  of  stored  air.  Demo  CRef.  6]  found  that  for 
the  supersonic  cascade  rig  in  building  230  an  area  of  6.27 
square  inches  was  needed  at  the  nozzle  throat  where  the  flow 
would  choke.   This  same  area  was  taken  for  the  total  outlet 
throat  area  of  the  swirl  vanes  for  the  high  speed  CDTD. 

Referring  to  the  following  sketch,  the  area  at  the  exit 
where  the  flow  will  choke,  A»,  is  given  by, 

A«  =  L   a   Z  <2) 

where  L  is  the  vane  length,  a  is  the  throat  width  and  Z  is 
the  number  of  vanes. 
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Figure  II. 1  Description  of  Conditions  at  Throat 

The  blade  spacing,  s,  is  given  by 

s  «  n  D  /  2  <3> 

where  D  is  the  diaseter  at  the  vane  exit.   Fros  equations 

<2>  and  (3), 

A*  -  L  (a/s)  n  D.  (4) 

In  Ref .  7,    the  following  conditions  were  calculated  for 

a  "full  scale"  high  speed  CDTD  design  to  have  a  test  section 

width  of  1  inch  at  25  inches  diameter: 

Mi»1.3     at  0i-55«,    at  Ri*12.5  inches 
Vi=1392  ft/sec    at  Tox520  °R,   giving 
V81-1140  ft/sec    and  Vr0*1357  ft/sec   and  for 
a  swirl  vane  length  of  5  inches,   Vp0=182  ft/sec 
and  B0a82.34°.   An  increase  of  swirl  vane  length 
to  5.713  inches  gave  Vro«160  ft/sec  and  Bo»83.3<» 

Using  the  above  data  for  the  "full  scale"  high  speed 
CDTD  in  equation  (4),  A«  «  38.5  inchea2.   Using  a  linear 
scaling  to  the  Model  in  which  A*  ■  6.27  inches?,  the  scale 
factor  is 
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(A-h/A*!)*  =  (6.27/38.5) *  %  .40. 


<3> 


Using  this  scale  factor,  the  diienalona  calculated  for  the 
high  speed  model  are  ehown  In  the  following  sketch. 
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Figure  II. 2  Croes  Sectional  View  of  the  Test  Section 

From  equations  (2), 

a  z  *  A«  /  L  (6) 

therefore  for  a  diameter  of  8  Inches  we  can  predict  a 
required  throat  dimension,  and  blade  spacing  using  equation 
(3)  for  a  selected  number  of  blades.   Results  are  given  In 
table  1.   The  requirement  then  Is  to  design  suitable  turning 
vanes  to  produce  the  required  throat  area  and  spacing. 
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III.  BLADING  DESIGN 

A.   BACKGROUND 

In  Ref.  5  Vavra  gave  a  procedure  to  establish  a  blade 
profile  for  specified  values  of  throat  width  a,  spacing  s, 
trailing  edge  thickness  te,  and  inlet  flow  angle,  <x.   Figure 
6  illustrates  these  features. 

Vavra' s  design  was  for  axial  turbine  blades  or  vanes.   A 
computer  program  was  available  to  construct  blading 
following  Vavra 's  methods  CRef.  5].   The  method  used 
circular  arcs,  and  straight  lines.   Wedge  angles  and  leading 
and  trailing  edge  radii  were  input  parameters  to  determine 
a   flow  channel  with  suitable  acceleration.   No  comparable 
method  and  program  were  found  for  radial  outflow  vanes. 

In  1966  Goldman  and  Scull  in  CRef.  83  designed  blades 
with  inlet  transition  arcs,  circular  arcs  and  outlet 
transition  arcs.   However,  these  blades  were  not  converging 
at  the  exit  and  were  not  designed  for  radial-outflow. 

Later  Goldman  and  Vanco  CRef.  93  developed  a  computer 
program  for  the  design  of  sharped-edged  throat  supersonic 
nozzles  for  isentropic  flow.   Their  program  was  constructed 
for  the  radial-inflow  case.   The  Goldman  and  Vanco  design 
also  involved  the  requirement  of  contour  machining. 

There  is  a  substantial  amount  of  information  on  both 
radial-inflow  and  axial  turbine  blade  design  in  the 
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literature   CRef.  8,9,10,11,  and  12  for  example].   Yet  none 
was  found  which  dealt  with  the  required  case  of  radial- 
outflow. 

B.   DESIGN  METHOD 

The  expansion  process  in  a  radial  turbine  differs 
appreciably  from  that  in  an  axial  turbine  because  of  the 
radius  changes.   Vavra's  work  CRef.  5]  involved  axial 
turbine  design.  The  work  presented  here  for  radial-outflow 
vane  design  followed  an  appreciation  of  Vavra's  methods  and 
possibly  retraced  steps  of  others  who  have  previously  sought 
the  simplest  geometries  for  machining  purposes. 

The  simplest  shape,  that  of  double  circular  arcs  was 
attempted  first.   Its  development  is  presented  in  Appendix 
B.   It  was  anticipated  that  with  the  requirements  of  radial 
inlet  flow  and  tangential  exit  outflow  the  double  circular 
arc  method  could  be  used  in  the  design  of  radial-outflow 
turning  vanes.   Once  the  method  was  developed  and  exercised 
it  was  clear  that  more  geometrical  constraints  were  required 
to  produce  a  nozzle  with  an  exit  which  was  both  nearly- 
tangential  and  convergent. 

The  double  circular  arc  program  generated  vanes  which 
appeared  to  converge  within  the  passage  and  then  diverge 
towards  the  exit  of  the  nozzle.   Therefore,  the  flow  would 
choke  within  the  nozzle  and  the  outlet  flow  conditions  could 
not  be  predicted  very  well. 
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Modifications  were  made  to  the  double  circular  arc 
design  to  achieve  convergency  at  the  exit.   The  first  step 
was  to  remove  the  pressure  line  and  use  the  camber  line  as 
the  pressure  surface.   The  assumption  was  made  that  the 
effects  associated  with  inlet  flow  incidence  resulting  from 
the  removal  of  the  pressure  line  would  be  negligible.   The 
results  indicated  that  the  passage  was  indeed  converging 
toward  the  exit.   The  second  step  in  an  attempt  to  assure 
tangency  and  convergency  at  the  exit  was  to  allow  for  a 
thickness  and  a  straight  section  at  the  leading  and  trailing 
edges.   The  straight  section  was   such  that  a  specified 
throat  dimension  was  achieved  at  the  trailing  edge  of  the 
adjacent  previous  vane.   The  resulting  program,  which  is 
presented  in  Appendix  C,  designs  a  wedge  shaped  arc  blade 
such  that  a  passage  which  converges  at  the  exit  is  realized 
with  a  specified  throat  dimension.   Using  the  method  a  high 
speed  CDTD  blading  design  was  obtained.   Figure  7  shows  the 
resulting  blade  design. 
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IV.  FLOW  ANALYSIS 

A.   METHOD 

A  FORTRAN  computer  program  for  solving  numerically  the 
two-dimensional  Euler  equations  using  the  Godunov  method  was 
developed  at  the  Turbopropulsion  Laboratory  by  Dr.  S. 
Eidelman  [Ref.  123.   With  Dr.  Eidelman's  help  the  program 
was  to  be  applied  to  obtain  a  prediction  of  the  flowfield 
through  the  vane  passage.   An  abbreviated  description  will 
be  given  here  of  the  computational  method. 

The  unsteady  two-dimensional  Euler  equations  are  written 
as 


where 
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^  =  density 
U  and  V  are  the  velocity  components  corresponding  to  the  x 

and  y  coordinates 
p  =  pressure 
T  =  specific  heat  ratio. 
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Equation  (7)  is  nondimensionalized  using  arbitrary 
reference  quantities,  and  the  Godunov  method  CRef .  12]  is 
applied  to  compute  at  each  time  interval  the  fluxes  at  the 
boundaries  of  the  cells  of  a  grid  established  for  the 
problem  under  consideration.   The  steady  state  solution  is 
found  when  no  further  change  occurs  between  successive  time 
steps. 

The  computational  domain  is  shown  for  a  typical  cascade 
problem  in  Figure  8  and  the  grid  and  cell  notation  are 
defined  in  Figure  9.   Referring  to  Figure  8,  for  the 
tangency  condition  at  the  wall  of  the  passage  a  solid  wall 
boundary  condition  is  imposed  over  the  segments  3-5  and  4- 
5.   Boundary  conditions  for  segments  1-3,  5-7,  2-4  and  6-8 
are  for  periodicity  of  the  flow  with  respect  to  the  y  axis. 

The  parameters  held  constant  at  the  inlet  are  the  x 
components  of  velocity,  U,  the  flow  angle,  9,  the  gas 
stagnation  enthalpy,  H,  and  the  entropy,  S.   The  velocity, 
density  and  pressure  can  be  calculated  from  these  quantities 
easily  as  V^n,  Pin*  anc*  Pin*   At  the  outlet  boundary  only 
pressure  is  defined  and  other  flow  parameters  are  calculated 
in  the  solution.   In  the  solution  procedure,  Uout,  vout>  and 
V  out  are  set  equal  to  the  values  of  U,  V  and  P  one  point 
ahead  of  the  outlet  boundary. 

A  subroutine  (program)  constructs  a  non-orthogonal  grid 
covering  the  computational  domain  as  shown  in  Figure  9.   The 
flow  parameters  related  to  the  center  of  the  cell  have  a 
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fractional  index:  j  +  h»  k*  5£.   The  cell  boundaries  have  one 
fractional  and  one  integer  index:  3,  k+ %   or  m+h,  n.   The 
parameters  at  the  time  t  will  have  subscript  indices  and  at 
time  t*  St  the  parameters  will  have  superscript  indices. 

The  Euler  equation  is  approximated  to  first  order  by  the 
equation  in  table  2.   Solution  to  the  Riemann  problem  is 
initiated  assuming  that  the  left  state  is  that   located  at 
3-%,k+h   and  the  right  state  is  that  at  j+ J$,k+ 5$. 

The  solution  to  the  Riemann  problem  will  give 

CR,P.Un.Vt,>3  ^ 

Un  =  the  velocity  component  normal  to 

the  edge  <j,k*l)  -<j,k) 

V*t  =  the  velocity  component  tangent  to 

the  edge  <j,k+l)  -<j»k) 

The  capital  letters  R  (rho),  U,  V,  P  are  the  parameters 

calculated  at  the  edges  of  the  zones. 

The  component  of  the  velocity  is  then  transformed  back 

to  the  cartesian  coordinate  system  and  a  final  value  of  the 

parameter  on  the  cell  edge  is  received 

<R,P,Un,B) 

3#k+  J* 

B.   RESULTS 

The  program  (EGE)  was  run  on  the  IBM  370-3033  at  the 
Naval  Postgraduate  School.   Figure  10  illustrates  the  grid 
covering  the  passage  for  the  wedge-arc  blading  shown  in 
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Figure  7.   An  enlargement  of  the  exit  area  is  shown  in 
Figure  11. 

In  order  to  apply  the  code  to  the  radial  geometry,  some 
modifications  were  necessary.   A  transformation  of 
coordinates  was  made  to  conveniently  accommodate  the  radial 
flow  inlet  boundary  condition.   Since  the  program  remained 
basically  two-dimensional,  the  outlet  flow  was  treated  as 
being  periodic  over  the  blade  space  at  the  exit  which  then 
did  not  change  with  radius.   Thus  the  solution  found  will  be 
valid  for  the  flow  within  the  rim  and  approximately  true 
just  outside  the  rim.   For  this  reason  the  computational 
domain  in  Figure  11  was  not  extended  radially  outwards 
beyond  one  blade  space. 

Results  could  not  be  obtained  for  the  flow  through  the 
passage  within  the  time  available.   Some  dif f iculities  were 
encountered  by  the  program  in  the  region  of  the  throat. 
Some  modifications  of  the  grid  must  be  made  before  a 
converged,  steady  flow  solution  can  be  obtained. 
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V.   SUMMARY  AND  DISCUSSION 

Two  exploratory  projects  at  the  Turbopropulsion 
Laboratory,  1).  to  teat  the  concept  of  detonative  combustion 
as  an  alternative  to  mechanical  compression  in  gas  turbines 
and  2) .  to  design  a  test  rig  for  radial  diff users  used  in 
centrifugal  compressors,  were  examined.   The  projects  were 
similar  in  nature  since  they  involved  the  generation  of 
nearly-tangential  radial  outflow  and  each  required  a  design  of 
turning  vanes  for  a  test  apparatus.   The  overriding  limitation 
for  either  project  was  that  any  apparatus  developed  would  have 
to  be  easily  and  inexpensively  manufactured  in-house. 

In  the  Rotary  Detonative  Wave  Engine,  "RDWE",  project, 
turning  vanes  were  required  such  that  upon  detonation  the 
maximum  possible  torque  impulse  would  be  obtained.   The  best 
design  possible  to  yield  maximum  power  hence  maximum  torque 
was  to  produce  turning  vanes  that  exit  nearly-tangential. 
Using  Euler's  Turbine  equation  a  relationship  was  developed 
relating  exit  angle  to  specific  power  and  efficiency. 

The  results  for  the  unsteady  expansion  process  showed 
that  the  value  of  the  relative  flow  angle  <B2>  established  the 
limiting  conversion  efficiency  and  specific  power.   If  the 
relative  velocity  increases  beyond  the  limiting  point  the 
specific  power  generated  will  be  absorbed  rather  than 
generated  and  efficiency  will  decrease.   High  efficiency 


29 


occurs  in  large  part  only  when  the  specific  power  is  low  or  at 
the  limiting  point.  In  general  the  conversion  efficiency  will 
vary  during  the  expansion  process,  but  will  always  be  poor. 

This  result  implies  that  the  concept  of  the  RDWE  can  not 
be  efficient  as  a  stand  alone  rotor.   Other  methods  would  have 
to  be  employed  to  increase  efficiency,  such  as  the  utilization 
of  additional  components,  particularly  those  which  would  allow 
nearly  constant  flow  by  sequencing  the  detonation  wave  with 
respect  to  partial  porta. 

The  second  project  involved  the  development  of  a  test  rig 
for  a  high  speed  Centrifugal  Diffuser  Test  Device,  "CDTD". 
The  project  would  closely  parallel  an  in-house  concept  of  a 
low  speed  CDTD.   The  low  speed  version  simulates  the  flow  of 
air  delivered  to  a  diffuser  as  a  steady  process.   The  high 
speed  version  simulates  in  a  steady  process,  the  transonic 
Mach  numbers  characteristic  of  advanced  centrifugal 
compressors . 

For  the  high  speed  CDTD  a  passage  would  be  required  which 
would  cause  the  flow  to  converge  and  choke  at  the  exit  of  the 
passage.   Because  the  concept  of  the  RDWE  was  questionable  the 
main  interest  tuned  to  designing  vane  shapes  and  the  apparatus 
for  a  high  speed  CDTD. 

Two  methods  and  their  respective  computer  programs  were 
developed.   The  first  method  involved  the  development  of 
double  circular  arc  vane  shapes.   Given  the  outer  radius,  the 
inner  radius  and  the  number  of  blades  the  radii  of  the  suction 
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and  pressure  sides  were  determined.   The  design  emphasis  was 
on  acheiving  a  tangent  suction  surface  at  the  exit  and  a 
radial  camber  line  at  the  inner  radius.   The  concurrent 
condition  of  a  tangent  camber  line  at  the  exit  gave  the 
solution  of  a  blade  of  zero  thickness. 

The  geometrical  description  resulting  in  a  tangent  camber 
line  at  the  exit,  resulted  in  the  sine  of  the  camber  angle  0 
being  equal  to  the  cosine  of  the  subtended  blade  angle  p    . 
Following  manipulations  to  relate  the  value  of  f*     to  just  the 
radius  ratio  (the  only  specified  parameter) ,  an  equation  of 
fourth  order  in  sin  r   resulted.   Using  Newton's  method  of 
interation  the  first  (smaller)  root  was  found.   Using  this 
value  of  i     a  blade  shape  was  generated.   It  was  clear  that 
the  suction  side  was  tangent  but  the  efflux  angle  a  was  not 
zero.   An  efflux  angle  of  zero  would  indicate  that  the  camber 
line  was  also  tangent  to  rim.    A  graph  of  the  fourth  order 
function  indicated  that  two  roots  were  possible  for  each 
radius  ratio.   Only  the  first  root  was  solved  for.   It  is 
probable  that  the  second  root  would  yield  the  tangency 
condition  of  both  camber  line  and  suction  side.   The  selected 
root  was  shown  however  to  give  sin  0  =  cos  f"   . 

It  was  realized  that  convergence  at  the  exit  could  not  be 
guaranteed  with  the  double  circular  arc  blade  design. 
Modifications  were  applied  to  the  design  to  cause  the  passage 
to  converge.   Improved  convergence  resulted  when  the  original 
camber  line  and  the  suction  line  were  used  to  describe  the 
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blade.   The  simple  double  circular  arc  design  required  curved 
or  thin  trailing  edges.   With  the  addition  of  a  finite 
thickness  at  the  leading  and  trailing  edge,  a  straight  line 
segment,  and  a  wedge  angle  all  of  which  would  better  direct 
the  flow  in  the  tangent  direction,  practical  blade  shapes  with 
convergent  passage  were  acheived. 

Geometrical  equations  were  developed  to  yield  the  proper 
blade  design.   The  exit  throat  dimension  was  fixed  in  the 
design  based  on  the  limitations  of  the  available  air  supply 
and  based  on  dimensions  scaled  from  a  high  speed  CDTD  design. 
It  was  anticipated  that  the  flow  would  choke  at  the  specified 
throat  dimension  between  adjacent  suction  and  pressure  sides 
and  the  straight  line  segment  would  intersect  the  exit  circle 
which  was  tangent  to  the  suction  side.   This  did  not  occur. 
The  length  of  the  straight  line  segment  varied  between  1  and  5 
percent  of  the  value  needed  for  this  condition.   It  was  not 
required  that  the  exit  circle  be  tangent  to  the  pressure  side. 

From  the  geometry  and  the  specified  parameters  a  blade 
was  developed  giving  convergency  at  the  exit.   The  camber 
angle  0  and  the  subtended  blade  angle  f*    were  specified  from 
experience  obtained  in  the  development  of  the  double  circular 
arc  design.   The  chosen  angles  0  and  j*    gave  an  efflux  angle  a 
between  8  and  9  degrees.   To  achieve  the  perfect  condition  of 
tangential  flow  angle  a  would  be  required  to  be  equal  to 
zero.   Based  on  current  technology  and  manufacturing  abilities 
an  angle  in  the  area  of  8  or  9  degrees  would  give  as  close  to 
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tangential  flow  as  is  possible.   It  is  noted  that  the  angle  0 
is  related  uniquely  to  the  angle  J*  •       Once  the  radius  ratio  is 
given. 

The  coordinates  of  the  designed  passage  were  then  input 
into  the  FORTRAN  program  (EGE)  to  predict  analytically  the 
flowfield  through  the  vane  passage,  by  solving  the  Euler 
equations.   The  program  (EGE)  was  designed  to  run  on  the  IBM 
370-3033.   The  program  successfully  constructed  a  non- 
orthogonal  grid  covering  the  domain,  however  results  could  not 
be  obtained  for  the  flowfield  through  the  passage  in  the 
initial  attempts.   It 

is  thought  that  the  difficulty  lies  in  the  curvature  and 
density  of  the  grid  in  the  region  of  the  exit  throat.   Since 
the  EGE  program  has  not  previously  been  applied  to  a  radial 
outflow  problem  with  severe  curvatures  the  task  of  analyzing 
the  present  flow  with  the  EGE  code  is  now  thought  to  require 
the  concentrated  effort  of  the  code's  author. 

In  contrast,  the  CDTD  vane  design  can  be  manufactured 
without  difficulty  and  the  flow  from  the  turning  vanes  can  be 
measured  experimentally.   It  is  suggested  that  such  an 
experiment  would  provide  a  severe  test  of  the  ability  of  the 
EGE  code  to  simulate  transonic  radial  outflow  containing 
embedded  waves. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

The  following  conclusions  are  offered: 

The  Rotary  Detonation  Wave  Energy  concept  in  its  original 
configuration  is  potentially  inefficient.   An  increase  in 
efficiency  would  require  the  addition  of  a  component  which 
could  recover  the  unsteady  kinetic  energy  of  the  flow  leaving 
the  rotor. 

A  program  for  designing  radial  outflow  vanes  using  double 
circular  arc  blades  was  demonstrated  successfully.   While  the 
method  may  be  useful  in  other  contexts,  it  was  found  to  be 
unable  to  produce  converging  passage  shapes  with  a  near- 
tangential  outlet  flow  direction. 

A  program  for  designing  radial  outflow  vanes  using  wedge- 
arc  blade  shapes  was  used  successfully  to  produce  converging 
passage  shapes  with  prescribed  throat  widths  and  with  a  near- 
tangential  outlet  flow  direction. 

Both  blade  shapes  can  be  produced  relatively  inexpensively 
from  standard  metal  stock. 

A  preliminary  design  of  a  test  model  of  a  high  speed 
Centrifugal  Diffuser  Test  Device  was  produced. 
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B.   RECOMMENDATIONS 

The  following  recommendations  and  areas  of  further 
investigations  are  suggested: 

An  arrangement  of  partial  porting  and  following  power 
turbine  for  the  Rotary  Detonation  Wave  Engine  should  be 
examined. 

The  program  for  double  circular  arc  blades  should  be 
modiied  to  seek  both  roots  of  the  condition  for  camber  line 
tangency,  and  to  offer  the  selection  of  any  value  of  Y>   between 
them.   This  will  result  in  the  ability  to  examine  any  value  of 
y^    within  the  range  of  any  practical  interest. 

In  the  case  of  wedge-arc  blades,  the  possibility  of 
solving  for  the  particular  case  in  which  the  circle  at  the 
passage  outlet  is  also  tangent  to  the  pressure  side  should  be 
examined.   A  procedure  should  also  be  introduced  to  calculate 
the  progression  of  the  passage  area  taking  into  account  the 
increasing  radius. 

The  proposed  blading  geometry  should  be  constructed  and 
the  outlet  flowfield  examined  experimently . 

Work  should  be  continued  to  obtain  a  computational 
prediction  of  the  flow  through  the  passage.   The  model  test 
could  then  serve  as  a  verification  of  the  code. 
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TABLE  1 .   Required  Design  Dimensions  for  Convergent 
Throat  in  High  Speed  CDTD. 


NUMBER  OF  BLADES         THROAT  DIMENSION  SPACING 

Z  a  (in.)  s  (in.) 

6  .5225  4.1888 

12  .2613  2.0944 

18  .1742  1.3963 

36  .08708  .6981 

45  .06967  .5585 

60  .05225  .4189 
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TABLE  2.   First  Order  Solution  of 
Euler  Equations 
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Outer  Plexiglas- 

Radius 


Teat  Section 


Note:   Not  to  scale. 


Figure  4.   CDTD  Flow  Generator  -  Section  AA 

from  Figure  3. 
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Figure  5.   Two-Dimonsional  Turbine  Ca 


ecada . 
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Figure  6.  Axial  Turbine  Blade  Geometrical  Design 
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Ro  *4    Ri  *  3-2    No.  Blades  »  60 

?>=    9.1.72521902   0  «  80.82738092 

Condition  of  Tangency  la 
sin  0  =  cos 
sin  0  =  0.159407755    cos    -  0.159409430 
chord  -  0.983540343   chord  tost  »  0.983540073 
alpha  »  18.34514098 
radial  solidity  =  1.909859317 
thickness  to  chord  ratio   -  0.194933797 
radius  of  caaber  circle    ■  0.758551242 
radius  of  suction  circle   -  0.575175292 
radius  of  pressure  circle  ■  1.462802226 

Figure  7.   Final  Blading  Design. 
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Figure  8.  Computational  Grid  and  Boundary  of  the 
Computational  Domain. 
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a>   Grid  Location  in  Coordinate  Syat 


e» 


K  +  l 


J    i*i    x 


b)   Single  Cell 


Figure  9.   Computational  Grid  Notation 
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Figur«  10.   Computational  Grid  for  Outflow 

Analysis 
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Figure  11. 


Enlargement,  of  the  Grid  in  the 
Throat  Region. 
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APPENDIX  A 
TORQUE  PRODUCED  BY  DISTRIBUTED  GAS  REACTION 

A.l   METHOD 

The  RDWE  concept  involves  an  unsteady  process.   The 
potential  performance  of  the  RDWE  can  be  understood, 
approximately,  by  examining  the  range  of  conditions  through 
which  the  unsteady  process  passes.   Thus  the  steady  flow 
performance  of  the  proposed  arrangement  needs  to  be 
examined.   We  first  attempt  to  predict  analytically  the 
power  or  torque  impulse  that  can  be  generated  from  the 
engine  under  steady  flow  outlet  conditions. 

The  passages  in  the  RDWE  exit  nearly  tangentially 
(Figure  1).   Since  this  condition  is  difficult  to  achieve  we 
examine  the  effect  of  the  exit  angle  on  the  amount  of  power 
that  can  be  generated  and  the  efficiency  that  can  be 
produced. 

The  basic  equation  for  relating  the  fluid  flow 
conditions  to  the  energy  transferred  to  the  shaft  is  Euler's 
Turbine  equation.  The  torque  of  the  fluid  acting  on  the 
rotor  is  given 

T  =  »  (  ri  V6i  -  T2  V02)  A<1) 

If  the  rotor  has  a  rotational  speed  w' ,  the  shaft  power 
transferred  or  rate  of  work  transferred  from  the  fluid  to 
the  rotor  is  given  by. 
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P   =  w'  T  A<2> 

The  work  per  unit  mass  flow  is  given  by 

P  /  m   =  w'  <ri  V9i  -  r2  V02)  A<3) 

The  subscript  1  denotes  entrance  station,  2  denotes  exit. 
Defining  the  rim  speed  as   U  =  rw',  equation  A<3)  becomes 

P  /  ra  =  Wx  =  Ui  Vei  -U2  V02  A<4) 

A  general  velocity  diagram  for  a  rotor  is  shown  in  Figure 
Al. 

The  tangential  component  of  velocity  V©  is  taken 
positive  in  the  direction  of  the  rim  velocity  U.   Figure  Al 
shows  all  positive  velocity  components  in  conventional 
cylindrical  coordinates. 

The  gas  enters  the  wheel  near  the  center,  effectively 
with  zero  moment  of  momentum  (about  the  axial  direction) . 
Therefore  V ei  =  0,  and  from  equation  A<4>, 

Wx  =   P   =   -U2  V02  A<5) 

m 

For  Wx  to  be  positive,  it  is  seen  that  V©2  *s  negative.   The 

exit  velocity  diagram  representing  the  flow  at  the  rim  of 

the  wheel  is  given  in  Figure  A2. 

From  the  geometry  shown  in  Figure  Al,  in  general 

sin  ft  =  V8  -  U  A<6) 

W 
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Figure  Al.   Velocity  Triangle  of  Conventional 
Cylindrical  Flow. 


H4— HfcjH 


Figure  A2.   Velocity  Diagram  at  Exit  of  Ri». 
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therefore 

V0  =  W  sin  B  +  U  A<7> 

Substituting  equations  A(7>  into  A(5)  we  obtain, 

Wx  =   P   =  -  U*2  C  W  sin  B  *    U  1 
m 

=   -  U22  C  W  sin  B  +  1  3 
U2 

U22  C  W  sin  <-B>  -  1  3  A<8> 

U2 

Defining  the  power  conversion  efficiency  "c  as  the 

fraction  of  the  available  power  actually  converted  to  shaft 

work,  thus 

»c  =     P/k A(9) 

P/m  +  V2  /2 

and  defining  B2'=  -B2»  the  relative  flow  angle  at  the  rim 
which  is  positive  in  the  direction  of  rotation,  then 

P   =  u2  [  W2  sin  B2'  -13.  A<10) 

m  U2 

From  the  geometry  shown  in  Figure  Al, 

V22  =  Ve22*  <  W*22  -CV02-U232) 
so  that, 

)4V22=  CU2  V62  -  W22+  W223  A  (11) 

Using  equations  A(5)  and  A(10) 

-  U2  V02  =  -Wx  =  U2  C_W2  sin  B2'  -  13  A<12) 

U2 

Substituting  equation  A(12)  into  A<11)  gives 
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JiV2    =     H   U22     C    -2    W?    sin    R2'     +    1    +       W?    2       3  A  (13) 

U2  U2 

Substituting  equation  A<10)  and  A<13)  into  A(9) 

ic   =      P/m 

P/m  ♦  HV22 

LU2  <W-/U->  ain  B_'  -  13 
_2   — 2 — "~2 2  

U2[<W2/U*2>sin  B2'  "  13  *  W22  C-2WsinR2'  ♦1  +   «2  23 


U2  U2 


=  (W?/U?>  sin  B?'  -  1 
-54  *     h    <W22/^22) 


»c  =  -  2  C  <W?/U?  sin  B2'  -  1  3  A(14) 

<  1  -  W2/U2) 


Equations  A(10)  and  A<14)  are  valid  for  a  discreet  jet, 
or  for  distributed  gas  at  the  rim.   For  the  distributed  gas, 
it  is  possible  to  calculate  the  referred  power  as  a  function 
of  wheel  radius,  giving  specific  power,  Pa.   Using  equation 
A<10), 

P3  =  P/mU2   =  [(W2/U2>  sin  B2'-  13  A<15) 

Therefore  equation  A (14)  becomes. 


2   P. 


<W2/U2>2  ~  1 

For  sin  B2'  =  1#  implying  that  B2'=90  degrees, 

'c        2 

<W2/U2>+1 
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Thus  the  conversion  efficiency  is  unity  when   W2/U2  =  1 
but  this  is  a  condition  which,  from  equation  A<15)  gives 
zero  power  output. 

A.  2   RESULTS 

The  value  of  the  specific  power  and  conversion 
efficiency  (given  by  equation  A(15)  and  A(16>>  as  a  function 
of  the  ratio  of  relative  flow  to  wheel  velocity,  for  various 
values  of  the  relative  exit  angle,  are  shown  in  Figure  A3 
and  Figure  A4  respectively. 

Figure  A3  shows  that  the  specific  power  increases 
linearly  as  the  exit  velocity  increases,  and  decreases  as 
the  relative  exit  angle  becomes  less  tangential .   The 
conversion  efficiency  is  seen  to  decrease  as  the  power 
increases  and  to  decrease  as  the  relative  flow  angle 
decreases  (less  tangential). 

The  results  shown  in  Figure  A3  and  Figure  A4  can  be 
used  to  gauge  the  performance  of  the  unsteady  RDWE.   In  the 
unsteady  expansion  process  the  value  of  the  relative  flow 
angle  (B2>  will  be  constant,  being  determined  by  the  passage 
exit  geometry.   The  wheel  speed  will  be  approximately 
constant  and  only  the  value  of  the  relative  velocity  will 
change.   Thus  the  unsteady  process  at  different  times  will 
be  described  by  points  on  a  single  line  shown  on  Figure  A3 
and  A4,  depending  on  the  particular  value  of  82* 
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The  difficulty,  clearly,  is  that  the  exit  velocity  must 
not  be  allowed  to  decrease  below  a  specific  factor  times  the 
wheel  speed  or  power  will  be  absorbed  rather  than  generated. 
Also  the  conversion  efficiency  will  vary  during  the 
expansion,  but  will  always  be  poor. 

This  implies  that  the  concept  of  the  RDWE  as  shown  in 
Figure  1  can  not  be  efficient.   The  device  would  require 
some  type  of  diffuser  that  would  convert  the  kinetic  energy 
leaving  the  rotor  into  static  pressure.   Unfortunately,  a 
variation  in  the  relative  velocity  magnitude  implies 
variation  in  both  magnitude  and  direction  of  the  exit 
absolute  velocity.   A  suitable  diffuser  for  such  a  flow  does 
not  exist.   Consequently,  other  geometries,  particularly 
those  which  would  allow  constant  flow  to  the  diffuser  by 
sequencing  the  detonation  wave  with  (partial)  ports  at  the 
rotor  exit,  should  examined. 
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TABLE  Al 


Specific  Power  and  Power  Conversion  Efficiency 


Vs 

Velocity 

Ratio 

Data 

60  DEG 

70  DEG 

w?/u? 

P« 

*c 

W-?/U? 

P« 

'c 

0 

-1 

2 

0 

-1 

2 

.5 

-.56699 

1.51200 

.5 

-.53015 

1.4137 

1.0 

-.13397 

oo 

1.0 

-.06031 

CD 

1.5 

.29904 

.47846 

1.5 

. 40953 

.65526 

2.0 

. 73205 

. 48803 

2.0 

.87938 

.58626 

2.5 

1.16506 

.44383 

2.5 

1.34923 

.51399 

3.5 

2.30109 

.36108 

3.5 

2.28893 

.40692 

4.5 

2.89712 

.30099 

4.5 

3.22862 

.33544 

5.5 

3.76314 

.25731 

5.5 

4.16831 

.28501 

6.5 

4.62917 

.22444 

6.5 

5.10800 

.24766 

7.5 

5.49519 

. 19892 

7.5 

6.04770 

.21892 

w?/u? 


80  DEG 
Ea 


*c 


90  DEG 


ic 


0 

-1 

2 

O 

-  1 

2 

.5 

-.50760 

1 

. 35359 

.5 

-.5 

1 

.33333 

1.0 

-.01519 

CD 

1.0 

0 

1 

1.5 

.47721 

.76354 

1.5 

.5 

.80000 

2.0 

.96961 

.64641 

2.0 

1.0 

.66667 

2.5 

1.46202 

.55695 

2.5 

1.5 

.57143 

3.5 

2.44683 

.43499 

3.5 

2.5 

.44445 

4.5 

3.43164 

. 35653 

4.5 

3.5 

.36364 

5.5 

4.41644 

.30198 

5.5 

4.5 

.30769 

6.5 

5.40125 

.26188 

6.5 

5.5 

.26667 

7.5 

6 . 38606 

.23117 

7.5 

6.5 

.23529 
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APPENDIX  B 
DOUBLE  CIRCULAR  ARC  BLADING  DESIGN 

B.l   METHOD 

A  computer  procedure  was  developed  to  design  radial 
out-flow,  near  tangential  turning  vanes,  specifying  only  the 
outer  radius,  the  inner  radius,  and  the  number  of  blades. 
The  blade  spacing,  radial  solidity,  camber  angle,  suction 
and  pressure  surface  radii,  chord  length,  thickness  to  chord 
ratio,  and  efflux  angle  are  determined.   The  suction  surface 
is  designed  to  be  tangential  at  the  rim. 

A  Hewlett  Packard  9830A  computer  was  used.   The  blade 
shapes  were  generated  using  a  Hewlett  Packard  9862A  plotter. 

The  geometrical  approach  to  the  design  of  the  vanes  is 
illustrated  in  Figure  Bl  and  much  of  the  notation  is  defined 
in  this  figure.   The  following  assumptions  were  made  as 
initial  constraints: 

1.  radius  at  exit  R0,  is  known 

2.  radius  at  entrance  Rj.,  is  known 

3.  number  of  blades  Z,  is  known 

4.  The  suction  side  is  tangent  to  the  outer 
radius  (point  C) . 

5.  The  camber  line  c',  is  tangential  at  the  inlet 
radius  (point  A) . 

6.  The  blade  subtends  an  angle  pat  the  center. 
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Initially,  a  procedure  was  developed  which  required 
specifying  the  camber  angle  0,  R0,  Ri,  in  addition  to  the 
number  of  blades.   Since  the  choice  of   0   which  would  give 
acceptable  values  of  efflux  angle  and  blade  thickness  was 
not  obvious,  an  equation  for  the  condition  that  the  camber 
line  was  tangent  at  the  rim  (which  would  give  a  blade  of 
zero  thickness,  since  the  suction  surface  was  also  required 
to  be  tangent)  was  written.   The  equation  was  seen  to  have 
two  possible  roots.   The  selection  of  the  smaller 
(incorrect)  root  gave  positive  values  of  the  efflux  angle 
and  reasonable  values  of  blade  thickness.   The  program  which 
was  developed  provides  the  option  of  selecting  0,  or  of 
selecting  the  blade  which  results  from  finding  the  incorrect 
smaller  root  to  the  tangent  camber  line  condition. 
Described  here  are  the  equations  for  the  generation  of  the 
blade  shapes. 

Defining  peripheral  blade  spacing  s  as, 

s  =  2  n  R0  /  Z 
and  the  angle  subtended  by  the  blade  space  as 

fa    =  360  /  Z 
the  radial  chord  as 

f  =  <R0  -  Ri>  /  R0 
then  the  radial  solidity  is  given  by 

<rr=(R0-Ri)/s=    P  /    ^a 
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Using  these  expressions  and  Figure  Bl  the  following 
equations  can  be  written: 

GE  =  Rc  -  Rc  cos  (0/2)  BCD 

=  Rc  CI  -  cos  (0/2) 
HE  =  c  T/2  +  Rc  C  1  -  cos  < 0/2]  B(2) 

Where  T  is  the  thickness  to  chord  ratio.  Also 

HE   =  Ra  C  1  -  cos  (0s/2)]  B<3) 

To  solve  for  the  thickness  and  T  we  subtract  the 
equations  B<1)  and  B<2) 

HE  -  GE  =  t/2  =  T  c/2  B(4) 

Using  Figure  Bl, 

FE   =  Rc  CI  -  cos  (0/2)3  -  T  c/2  B(5) 

=  Rp  CI  -  cos  (0p/2)3  B(6) 

Subtracting  equation  B(5)  from  B(l), 

GE  -  FE  =  t/2  =  T  c/2  B(7) 

The  chord  c,  is  given  by, 

c  =  2  Rp  sin  (0p/2)  B(8) 

=  2  Rs  sin  ( 0s/2)  B(9) 

=  2  Rc  sin  ( 0  /2)  B(10) 

All  equations  were  made  dimensionless  by  dividing  by  R0, 
giving  for  example 

c/R0  =  c"  B(ll) 

C"  =  2  Rp.  sin  (  0p/2)  B(12) 

Ro 

=  2  Ra  sin  ( 0s/2)  B(13) 

Ro 
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=  2  Re  sin  <  0  /2)  B<14) 

Ro 

For  the  geometry  we  seek,  given  RQ#  and  the  camber 
angle  0  and  efflux  angle  a,  the  angle  ACO  (or  ft)  is  given  by 

ft  =  90  -  (  a  +  0/2)  B<15) 

From  the  triangle  ACO 

f  *    90  -  (a  ♦  0/2)  ♦  180  -  0/2  =  180 

resulting  in, 

90    *  ^       -     0    *     Oi  B(16) 

This  equation  fully  defines  a.   The  values  of  ^  and  0 
are  needed  to  complete  the  description  of  the  camber  line. 
To  solve  for  the  radius  of  each  blade  arc  the  following 
relations  were  used. 
First  with 


As   =  Ec    -    *bz    cos  0/2  =  HE  B<17) 

Ro    Ro 

and 


Ap  =  Rp  -  Rp  cos  <0p/2)  =  FE  B<18) 

Ro    Ro 

Using  equation  B<3),  B<6)  and  B(7) 

Rc  C  1  -  cos  (0/2)3  -  Rp  C  1-cos  <0p/2)3  =  T  c/2    B<19) 

Also, 

A3  -  T  c/2  -  Rp  CI  -  cos  (0p/2)3  =  T  c/2         B<20) 
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which  gives 

Ap  =  As  -  T  c  B<21) 

Using  this  relationship  between  As  and  Ap  we  can  solve  for 
Rs  and  Rp.   From  equation  B<2>,  and  B<3>, 

As  =  Rs  CI  -  cos  <0s/2)  B<22) 

or 

cos  <0s/2)  =  1  -  As/Rs 
Squaring  and  adding  together  equations  B(22),  and  B(10) 
we  obtain 

cos2  <0s/2)  ♦  sin2  < 0s/2>  =  (1  -  As/Rs>2  +  <c/2Rs)2 
or 

1=1-  2<AS/Ra)  ♦  <AS/R3)2  ♦  C2/4RS2 
which  becomes 

Rs  =  J$AS  ♦   c  2_1 B<23) 

8Aa 

or 

Rfl  =  hAs  *   c  2  1 B<24) 

Ro  ^o   BAg 

Using  equations  B<5),  B(6),  and  B<16)  a  similar  manipulation 
will  yield, 

Rp  =  Mp  +  c2   1  B(25) 

Ro  SAp 

The  angles  0s  and  0p  are  found  using  equations  B(13), 

and  B(14)  since 

sin  <0s/2)  =  c/2Rs 
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or 


0s  =  2  arc  sin  <  c/Ro)  B<26) 

2RS/R0 


and  similarly 

sin  <0p/2)  =  c/2Rp 

or 

0p  =  2  arc  sin  <  c/Rn)  B<27) 

2Rp/RQ 

The  desired  condition  is  that  the  suction  surface  be 

tangent  to  the  outer  rim.   From  Figure  Bl  this  is  given  by 

the  relation 

0s/2  =  a  +  0/2  B<28) 

T  is  not  an  independent  parameter  since  it  can  be 

calculated  using  equation  B(25),  thus 

T  =  C  (Rg  C  1  -  cos  (0S/2)]  ■  Rr  I  1  -  cos  (0/2)]   )    B<29) 

<c/2) 

To  insure  near  tangential  flow  at  the  exit  we  first 

seek  the  condition  of  a  tangent  camber  line  at  the  exit. 

This  implies  since  we  also  seek  a  tangent  suction  surface, 

that  the  blade  would  have  zero  thickness.   Using  the 

notation  in  Figure  B2  the  following  equations  apply: 

sin  (  0/2)  =  L  /  c  BOO) 

and 

sin  f>    =  L/  R0  B(31) 

Substituting  equation  BOO)  into  BOD  gives, 

c  =  Rr,  sin  1£  B<32) 

sin  (0/2) 
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From  the  geometry 

cos  ( 0/2)  =  y/c  B(33) 

and 

cos  f      =  Ri     *    v  B(34) 

Ro 

or 

Ro  cos  f      =  Ri  +  y  B(35> 

Substituting  equation  B<33)  into  B<35)  gives, 

R0  cos  f      =  Ri  ♦  c  cos  <0/2)  B<36) 

or 

Ro  cos  'P      -  Ri  =  c  cos  (0/2)  B(37) 

Equation  B<32)  can  be  rewritten  such  that, 

R0  sin  ¥     =  c  sin  (0/2)  B<38) 

and  dividing  equation  B<38)  by  B(39), 

Ro  sin  jf     =   tan  (0/2)  B<39) 

Ro  cos  f      -    Ri 

This  equation  establishes  the  relationship  between  ^ 

and  0.   Using  Figure  B3  the  relationship  of  c,  r  ,  0  is 

established, 

cos  (0/2)  =   c/2  B<40) 

R  tan  f> 

or 

c  =  2  Ro  tan  f     cos  (0/2)  B(41) 

This  is  the  condition  for  the  camber  line  to  be  tangent  at 

the  exit. 

Using  equations  B<32)  and  B(41), 
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sin  P    =  2  tan  £•  cos  (0/2)  B(42) 

sin<  0/2) 

giving, 

cos  ^  =  2  sin  (0/2)  cos  (0/2) 

=  2  C  J*  sin  <  0/2  ♦  0/2)  ♦  sin  <  0/2  -  0/2)  ] 

■  sin  0  B(43) 

Therefore  for  the  camber  line  tangency  condition 

cos  f    =  sin  0  B<44) 

Using  equation  BOS)  and  B(44)  we  obtain  the  following 

2  sin  <  0/2)  cos  ( 0/2)   =  tan  <  0/2) 
2  cos  2  <  0/2) 

and 

2  sin  (0/2)  cos  (0/2)   =    Ro  sin  Z_ 

2  cos  2  (  0/2)         Ro  cos  y-      -    Ri 

Using  the  general  trignometrical  relationship, 

sin  2x  =  2  sin  x  cos  x 

we  can  manipulate  the  preceeding  equation  to  obtain 

sin  0  C  Ro  cos  ^   -  R<  ]  =  2  cos  2  (0/2)       B<45) 
Ro  sin  -J* 

Using  equation  B(43)  this  becomes 

cos   ^   C  R0  cos  -V  -  Ri  1       =2  cos2  (0/2)      B(46) 

Using  the  relationship 

cos2  x  =  sin2  x  +  cos  2x 

2  cos2  (0/2)  =  2  Csin2<0/2)  *  cos  0] 

=  2  C  1  •*•  cos  0  1 
2 
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=  1  +  (1-  sin2  0)  rL 
Using  equation  B(43)  this  becomes, 

2  cos2  (0/2)  =  1  ♦  (1  -  COS2  0)*4 

=  1  ♦  sin  ^  B<47) 

This  gives 

R0  cos2  y   -  Ri  cos  Y      =  R0  sin  ^   +  RQ  sin2^        B<48) 

Writing  s  =  sinO^/  equation  (48)  becomes, 

R0<1  -  s2)  -  Rj.  (  l-s2)'5  =  R0  s  >  R0  s2  B(49) 

or 

Ro  C  <l-s2)  -  s2  -  a]  =  -(1-S2)3* 
Ri 

(Rn)2   [  l  -  2s  -  3s2  ♦  4s3  +  4s<*  ]  =  l-s2 
<Ri> 

Thus  equation  B<48)  gives  the  condition 

F<s)  =  1Rj_>2  <l-32)  -  1  ♦  «(2  ♦  3a  -  4s2  -  4s3  )  =  0    B<50> 
Ro 

The  condition  of  tangency  of  the  camber  line  is  given 
by  the  solution  of  equation  BOO)  for  s  between  0  and  1, 
since  s=sin  /r    .   Using  Newton's  method  of  interation  the 
root  can  be  found  for  a  given  radius  ratio.   Figure  B4  shows 
F(s)  in  equation  B(50)  for  various  radius  ratios.   Figure  B4 
clearly  identifies  two  possible  roots.   The  larger  root  is 
the  correct  root,  giving  zero  blade  thickness.   The  smaller 
root  if  incorrect,  but  gives  positive  values  of  a  and  finite 
blade  thickness.   Emphasis  was  placed  on  finding  the  first 
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root.   The  value  found  will  give  the  selected  value  of  sin^. 
Once  ¥"     is  found,  using  equation  B(39)  0  can  be  calculated 
and  c  is  obtained  using  equation  B<32). 

Using  equation  B<28)  and  B<9)  0s  and  Rs  can  be 
calculated  and  the  thickness  is  then  given  using  equation 
B<29> . 

A  test  is  made  to  assure  that  the  length  Ap  is 
positive.   If  Ap  is  negative  this  implies  a  convex  rather 
than  concave  pressure  side.   This  would  also  imply  that  the 
inner  radius  is  too  small  and  should  be  increased  to  produce 
concave  curvature  on  the  pressure  side.   Ap  is  found,  from 
equations  B<19)  and  B(16). 

The  alternate  procedure  did  not  require  the  camber  line 
to  be  tangent  at  the  exit.   In  this  case  the  value  of  y   is 
not  solved  for  but  is  left  as  an  input  parameter.   The 
values  of  0  and  c  are  then  found  directly  using  equation 
B<39) . 

The  centers  of  curvature  are  given  by. 


*ds  =_Ec  _  <_R_c  cos  (0/2)  -  Rfl  cos  < 0s/2)  )  cos  (0/2) 

Ro  ^o              ^o 

Yds  =_B_L  ♦  <_Bc  cos  (0/2)  -  Ro  cos  (0s/2)  )  sin  (0/2) 

Ro  ^o              ^o 

xdP  =  Ec  ♦  <  Rp  cos  <  0p/2)  -  Re  cos  (0/2)  )  cos  (0/2) 

Ro  ^o               ^o 

ydP  =  B_L  "  <  EP_  cos  (  0p/2)  -  Re  cos  (0/2)  )  sin  (0/2) 

Ro  ^o               Ro 


with 
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Xd  =  Rc 

=  R^/Ro 

yd  =  Ri 

=  Ri/R0 
To  draw  successive  blades  on  a  rotor  the  program  must 
rotate  the  axes  about  the  rotor  center  after  each  blade  is 
drawn.   The  following  equations  were  used  (see  Figure  B5) . 

X  =    _x ♦  Y  sin  9 

cos  9      cos  © 

Y    ♦   x  sin  9  =  y 


cos  9       cos  9 
therefore, 

X  =  x  cos  9  +  y  sin  9 

Y  =  y  cos  9  -  x  sin  9 
Given  the  number  of  blades,  the  following  operations 
(for  all  three  arcs)  are  performed,  the  coordinates  are 
calculated  as 

x  =  xd,ds»dP  -  R  coa  ® 

y  ■  Yd>ds>ds  +  R  sin  ® 
The  arcs  are  then  drawn  from  9  =  (  (0/2)  -  (0s,p/2)  )  to 
((0/2)  +  (0s,p/2)) 
using 

X  =  x  cos  (N)  360  ♦  y  sin  (N)  360 

Z  Z 

and 

Y  =  y  cos  (N)  360  -  x  sin  (N)  360 

Z  Z 

where  N  is  the  number  of  blades.   Figure  B6  and  B7 
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illustrates  the  computer  generated  shapes  and  passage  for  a 
60  blade  design. 

B.2  COMPUTER  PROGRAM  "ACIRC" 

B.2.1   Description 

A  computer  program,  "ACIRC"  was  written  to 
construct  double  circular  arc  blading  using  the  above 
procedure.   The  program  will  calculate  and  draw  blading  for 
a  initially  radial-outflow  using  Newtons  Method  to  solve  for 
a  suitable  value  of  i*  ,  or  for  an  arbitrary  input  value  of  f  . 
In  either  case  the  method  will  yield  a  suction  surface  which 
is  tangent  to  the  outer  rim. 

The  program  is  user  friendly  and  is  annotated  for 
ease  of  understanding.   The  program  is  divided  into  six 
major  parts. 

The  first  part  is  the  input  data  section.   The 
second  part  calculates  the  periodicity,  radial  solidity, 
spacing,  and  asks  the  user  if  the  value  of  ^  is  to  be 
calculated  or  to  be  specified. 

The  third  part  calculates  the  first  root  for 
tangency  conditions  at  the  exit.   No  input  is  required  for 
this  section.   The  fourth  part  calculates  the  parameters 
needed  for  the  suction  surface  to  be  tangent  to  the  outer 
rim.   Usually,  if  the  user  specified  a  value  of  ¥      there  is 
the  possiblity  of  the  distance  from  the  chord  to  the 
pressure  circle  being  negative.   If  this  occurs  the  result 
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is  a  convex  arc.   To  remedy  this,  the  program  asks  the  user 
to  increase  the  dimension  on  the  inner  radius  by  a 
percentage.   The  program  will  use  this  value  to  solve  for 
the  required  inner  radius  to  achieve  a  suction  surface  that 
is  tangent  to  the  outer  rim. 

The  fifth  part  prints  the  output.   The  sixth  part 
draws  the  blade  surfaces. 

B.2.2   Required  Inputs 

The  following  inputs  are  requested: 

Outer  radius 

Inner  radius 

Total  number  of  Blades 

Yes  or  No  response  for  solving  for  the  value  of ^ 

Percentage  increase  in  the  inner  radius  -  may 

be  required,  when  the  value  of  ^  is  specified 
Number  of  Blades  requested  to  be  drawn 

B.2.3   Listed  Outputs 

The  following  outputs  are  listed: 

Outer  radius 

Inner  radius 

Number  of  Blades 

Angle  y» 

Angle  0 

Results  of  test  for  tangency 


71 


Chord  length 

Angle  a 

Thickness  to  chord  ratio 

Radius  of  camber  line 

Radius  of  suction  line 

Radius  of  pressure  line 
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B.2.4   Prograii  Listing 

Listed  hare  ia  the  program  listing, "ACIRC"  for 
the  design  of  double  circular  arc  blades 

10  REN*******ACIRC  1  JUNE  1934********** 

20  REM*******PR0GRRN  TO  CALCULATE  $<  DRAW  ARC  BLADING******* 

39  REM******F0R  AN  INITIALLY  OUTFLOW  CASCADE****** 

48  REN******U3ING  NEWTON'S  METHOD  TO  SOLVE  FOR  THE******** 

50  REM****0F  PSI********* 

60  DEG 

70  SCALE  -i» li-l» 1 

SO  DISP  "SET  PLOTTER  7  00  BY  7< Y>-<C0NT> "  5 

90  STOP 

95  REM*****INPUT  DATA************************ 

100  DISP  "ENTER  THE  OUTER  RADIUS"; 

HO  INPUT  RO 

120  DISP  "ENTER  THE  INNER  RADIUS"? 

ISO  INPUT  Rl 

140  R4=R1/R0 

150  DISP  "ENTER  THE  NUMBER  OF  BLADES"; 

160  INPUT  ZO 

170  L=0 

189  REM********END  INPUT  DATA**** 

1S5  REM*****CALCULATIONS********************* 

190  S5  =  360--'Z0 

200  PI =3. 141592654 

210  S2=2*P1*R0/Z0 

215  R9=(R0-R1)/S2 

220  S1=(R0-R1)/S2 

230  S0=S2-'R0 

240  R2=S0*S1 

25©  DISP  "WANT  PSI  CALCULATED  Y=1»0=N"5 

260  WAIT  10O0 

270  INPUT  01 

2S0  IF  01=1  THEN  340 

290  DISP  "INPUT  PSI"; 

300  INPUT  S3 

310  02=C0SS3 

320  T0=ATN(O2/SQR(1-O2t2)> 

330  GOTO  330 

340  GOSUB  1350 

350  REM*****N0W  DETERMINE  THE  VALUES  FOR  PSI, PHI? AND  C***** 

360  S3=ATN(K9/SQR(1-K9T2>) 

370  T0  =  2*ATN(R0*SINS3/''sR0*COSS3-Rl)) 

380  05=SINS3 

390  O6=COST0 

4O0  C1=R0*SINS3/SIN<T0'2) 

405  REM************************************** 

410  PEM******TEST  FOR  TANGENCY***** 
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420 

436 

440 

450 

460 

470 

4S0 

435 

490 

500 

510 

520 

530 

540 

550 

560 

570 

530 

590 

600 

6 1 0 

620 

630 

6  4  0 

650 

6  6  0 

670 

630 

690 

695 

700 

710 

720 

730 

740 

750 

760 

770 

730 

790 

868 

310 

320 

830 

840 

845 

850 

866 

876 

830 

896 

900 


>-<R3*<  l-C03<T0/2> )  > 
:0*Tl/2 


IS  TOO  SMALL" 5 
I.E. =25*  ENTER  H"5 


C2=2*TflNS3*COS<T0/2)*R0 

REM******CONCLUDES  CALCULATION  FOR  PSI » PHI i AND  C******* 

A0=90+S3-T0 

C0=C1/R0 

R3=C0/(2*SIN(T0/2>) 

R4=R1/R0 

D0=ATN<R3/R4> 

REM****FQR  SUCTION  SURFACE  TANGENT  TO  OUTER  RIM**** 
T2=<2*A0)+T0 

R5=C0/<2*SIN<T2/2>> 

T  1  =  ( <  R5*  ( 1  -COS  C T2 •■•"2 )  >  >  - <  R3*  < 1 -COS  ( TO  ••'2  > )  > )  ■•■"  <  CO .-'2 ) 

R1=R3*C1- C 0 S < T 0 ••••' 2 ) ')  + 1 

A2=A1-<C0*T1  > 

IF  A2>0  THEN  640 

IF  L=l  THEN  620 

DISP    "NO    SOLUTION-RI 

IUSP    "  "5 

DISP  "H  •  INCREASES  RI 

INPUT  H 

L=l 

R1=R1+(H/100) 

GOTO  200 

R6=(A2/2>  +  (  C0t2/>:S*A2>> 

P4=C0/(2*RS) 

T3=2*ATN(P4/SQR(1-P4t2>) 

0  3  =  R  3  *  R  0 

0  4 = R  5  *  R  0 

O1=R6*R0 

F  E  M  ^*^^  +  ^*^  +  ******^**^*  +  +  +  *^**^*****'^'A-","'i,r 

REM*******PRINT  STATEMENTS*************** 

PRINT  "P0="RG"RI="R1"N0.  BLADES="Z0 

PRINT  "PSI="S3 

PRINT  "PHI="TO 

PRINT  "CONDITION  FOR  TANGENCY  IS" 

PRINT  "     SIN  PHI=COS  PSI" 

PRINT  "   SIN  PHI  ="05"   COS  PSI  ="06 

PRINT  "CHORIi  =  " CI" CHORD  TEST="C2 

PRINT  "ALPHA="AO 

PRINT  " RADIAL  SOLIDITY=" > R9 

PRINT  "THICKNESS  TO  CHORD  RATI0="»T1 

PRINT  "RADIUS  OF  CAMBER  CIRCLE  ="03 

PRINT  "RADIUS  OF  SUCTION  CIRCLE  ="04 

PRINT  "RADIUS  OF  PRESSURE  CIRCLE ="01 

REM******END  PRINT  STATEMENTS***** 

F' EM  ************************************** 

REM****CENTERS  CAN  NOW  BE  DETERMINED**** 


Xl=R3-((R3*C0S(T0/2>)-(R5*C 
Y 1  =R4  +  (.  ( R3*C0S < TO/2  >  )- < R5*C 
X2=R3+ ( ( R6*C0S <  T3/2  > ) - ( R3*C 
Y2  =  R4-<<R6*C03(T3--2)  )-(R3*C 
X3=R3 


0S(T2/2)>  >*C0S(T0-'2> 
OS ( T2/2 >  >  )  *S  I N ( T0-"2 > 
OS ( TO/2  > ) >  *COS (T0/2 ) 
OS<T0/2>>  >*SIN<T0/2> 
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910 
929 
930 
940 
950 
9  6  0 
970 
9  8  0 
9  9  0 

0~ 


10 

20 

30 

40 

050 


Y3=R4 

REM****HRAW  BLADE  SURFACES**** 

IUSP  "INPUT  THE  NUMBER  OF  BLADES  TO 

WRIT  1000 

DISP  "      "5 

DISP  "STARTING  AT  N=l"5 

INPUT  N 

REM EDIT  NEXT  LINE 

N9  =  25 

FOR  N0=1  TO  N 

I  1  =  '::T0/2>-(T2/2> 

1 2  =  <  T  8  ••■''  2  >  +  (  T  2  ■  ■'  2  > 

I  3= (I 2- I 1>  N9 

REM***FOR  THE  SUCTION  SIDE*- 


DRAW 


FOR  PLOT  RESOLUTION' 


FOR  T4=I1  TO  12  STEF 
060  X4=X 1 - ( R5*C0S (T4) ) 


13 


•370 
6S0 
0  9  0 
100 
110 
120 
130 
140 


r4=Yl+<R5*SIN<T4)> 

X9=X4*C0S ( ( NO- 1 ) *S5> + ( Y4*SI N ( ( NO- 1 > 

Y9=Y4*C0S  < < NO- 1 >*S5 > - <  X4*S I N  < < NO- 1 ) *S5 ) > 


S5)) 


Y9 


T4 


PLOT 
NEXT 
PEN 

REM****FQR  CAMBER  LINE**** 

FOR  T4  =  0  TO  TO  STEP  T0.-N9 
1 5  y  X  4  =  X  3  -  '•  R  3  *  C  0  S  < T  4 ) ) 
160  Y4=Y3+CR3*SIN<T4>> 

1 7  0  X  9  =  X  4  *  C  0  S  <  (.  N  O  - 1  >  *  S  5  >  +  <  Y  4  *  S I N  <  < N  0  - 1  >  *  S  5  >  > 
1  SO  Y9=Y4*C0S ( ( NO- 1 ) *S5  > - ( X4*S I N ( ( NO- 1 ) *S5) > 
190  PLOT  X9jY9 
200  NEXT  T4 
210  PEN 

220  REM**#*FOR  PRESSURE  SIDE**** 
230  I4=T0/2-T3/2 
240  I5=T0/2+T3/2 
250  I6=(I5-I4)/N9 
260  FOR  T4=I4  TO  15  STEP  16 
270  X4  =  X2-(R6*C0S(T4) ) 
2S0  Y4=Y2+(Rb*SIN(T4>> 

290  X  9  =  X  4  *  C  0  S  ( < N  0  -  1 )  *  S  5 ')  +  Y  4  *  S I N  ( ( N  0  - 1 )  *  8  5  > 
300  Y9  =  Y4*C0S'::  (NO-1  )*S5>-X4*SIN(  (NO-1  >*S5> 
305  PLOT  X9>Y9 
310  NEXT  T4 
320  PEN 
330  NEXT  NO 

PEM**********"'"**",'**********':'"****",""'":':'**'"*-** 

END 

REM*****NEMTONS  WILL  BE  USED  TO  CALCULATE**- 

PEM*****TANGENCY  ROOT  BETWEEN  0  &  90  DEGREE' 

PEM*****SOLVING  FOR  THE  1ST  ROOT********** 

T9=0. 00001 


■->  O  -J 

340 
350 
360 
370 
3S0 


;**■*■** 
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1390  K9=@. 1 

1 400  DEF  FHF< K9 )  =  < 1-R4*2 > * < 1 -K9t2) -K9* <2+2*K9-4*K 

1410  DEF  FND<K9>=-<6*K9-2*K9*R4t2+12*K9t2+16*K9t3 

1420  K1=K9 

1430  K2=FNF<K9)/FND<K9> 

1440  K9=K1-K2 

1450  IF  flBS<K2)>fiBS<K9*T9)  THEN  1429 

1460  RETURN 


9t£-4*K9t3) 
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ACIRC" 

te: 

AO 

oc 

Al 

As 

A2 

Ap 

CO 

c" 

TABLE  Bl .   List  of  Computer  and  Text  Symbols 
for  Double  Circular  Arc  Blades 


IDENTIFICATION 

efflux  angle 

suction  line  to  chord  length 
pressure  line  to  chord  length 
c/R0  dimensionless  parameter 
for  chord  length 

CI         c  chord  length 

C2         c  chord  length  for  tangent  test 

11  -  Interval  ( 0/2  -  0s/2) 

12  -  Interval  (0/2  ♦  0s/2) 

13  -  Interval  (0/2  +  0s/2)  -  (0/2 

-  0s/2) 

Kl         -  constant  value 

K2         -  Newtons  method,  y/y' 

K9         -  initial  guess  for  Newton's 

method 

L         -  counter 

N         N  Number  of  blades  to  draw 

NO         -  Number  of  blades  to  draw  for 

loop 

N9         -  Number  of  points  for  interval 

of  arc 

01         -  Equals  1  for  solution  of  f  , 2 

for  input  "J£ 
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TABLE  Bl.   List  of  Computer  and  Text  Symbols  for 
Double  Circular  Arc  Blades  (con't) 


'ACIRC 


TEXT 


IDENTIFICATION 


02 

cos  ^ 

03 

Re 

04 

Rs 

05 

sin 

06 

COS   0 

PI 

n 

P4 

- 

R0 

Ro 

Rl 

Ri 

R2 

R3 

Rc^Ro 

R4 

Ri/Ro 

R5 

Rg/Ro 

R6 

Rp/Ro 

SO 

ys 

SI 

*r 

S2 

s 

S3 

t 

S5 

i* 

TO 

0 

cos  y- 

output  print  of  Rc 

output  print  of  Rs 

for  tangent  test 

for  tangent  test 

n 

constant,  cRo/2Rp 

outer  radius 

inner  radius 

density,  (R0  -  Ri>/R0 

camber  radius,  d  intension  less 

inner  radius,  dimensionless 

suction  radius, dimensionless 

pressure  radius, dimensionless 

periodicity,  2nRo/Z 

radial  solidity 

spacing 

blade  curvature  with  respect 
to  camber  line 

periodicity,  360/Z 

camber  angle 
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TABLE  Bl.   List  of  Computer  and  Text  Symbols  for 
Double  Circular  Arc  Blades  (con't) 

IDENTIFICATION 


thickness  to  chord  ratio 

suction  angle 

pressure  angle 

loop  parameter  to  draw  arcs 

tolerance  on  Newtons  method 
of  integration 

suction  center  x  coordinate 

pressure  center  x  coordinate 

camber  center  x  coordinate 

variable  for  x  values  in 
continuous  loop 

X9         -  variable  for  x  values  in 

rotational  loop 

suction  center  y  coordinate 

pressure  center  y  coordinate 

camber  center  y  coordinate 

variable  for  y  values  in 
continuous  loop 

Y9         -  variables  for  y  values  in 

rotational  loop 

ZO         Z  number  of  blades 


•ACIRC" 

TEX 

Tl 

T 

T2 

0S 

T3 

0p 

T4 

e 

T9 

- 

XI 

xds 

X2 

*dP 

X3 

*d 

X4 

— 

Yl 

Yds 

Y2 

YdP 

Y3 

Yd 

Y4 

_ 

79 


Figure  Bl.  Double  Circular  Arc  Blad< 
Geometrical  Description 
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•R,tu^ 


Figure  B2 


Geometrical  Description  for  Camber  Line 
Tangency  Condition. 
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Figure  B3.   Continued  Gooaatrlcal  Description  for  Camber 

Line  Tangency  Condition 
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Figure  B4.   Function  of  y7  Resulting  froi  the  Condition 
of  Tangency  of  the  Camber  Line  (equation  B<50>). 
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Figure  B5.  Schematic  for  Rotation  of  Axes. 
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Ro  =4    Ri  "  3-2    No-  BlAdaa  ■  60 

y*  9.1.72521902   0  =  80.82738092 

Condition  of  Tangency  is 
ain  0  =  coa 
ain  0  =  0.159407755    coa    *  0.159409430 
chord  =  0.983540343   chord  teat  »  0.983540073 
alpha  »  18.34514098 
radial  aolidity  =  1.909859317 
thickneaa  to  chord  ratio   *  0.194933797 
radiua  of  canber  circle    ■  0.758551242 
radius  of  auction  circle   «  0.575175292 
radiua  of  pressure  circle  «  1.462802226 


Figure  B6.  Computer  Generated  Shapei 
for  a  60  Blade  Design. 
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APPENDIX  C 
WEDGE-ARC  BLADING  DESIGN 

C.l   METHOD 

Because  of  the  limitations  seen  in  the  outlet  passage 
shapes  which  could  be  generated  using  purely  double  circular 
arc  blades,  symmetrical  linear  surface  segments  were  added 
in  an  attempt  to  achieve  a  convergent  nozzle  at  a  shallow 
angle  at  the  rim  exit. 

From  the  design  using  only  double  circular  arcs  it  was 
found  that  by  using  the  original  camber  line  and  the  suction 
line  to  describe  the  blade,  improved  convergency  at  the  exit 
was  obtained  (see  figure  CD.   Removing  the  original 
pressure  line  would  cause  the  flow  to  be  at  a  small 
incidence  angle  at  the  inlet,  but  this  was  thought  to  be 
acceptable. 

By  allowing  the  flow  to  exit  over  a  straight  suction 
side  from  the  exit  of  the  passage,  better  control  of  the 
outlet  flow  angle  was  assured.   The  wedge  arc  geometry  is 
shown  in  figure  C2. 

Referring  to  figure  C2,  the  wedge  arc  blade  shape  was 
developed  to  have  the  following  features,  or  constraints: 

1.  The  pressure  side  is  a  circular  arc. 

2.  The  suction  side  between  point  H  and  I  is  a  circular 
arc,  between  straight  and  symmetrical  line  segments. 
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3.  A  wedge  angle  <()  is  specified  at  the  leading 
and  trailing  edges. 

4.  Equal  leading  edge  and  trailing  edge  thicknesse  ( S) 
are  specified. 

5.  The  radius  to  the  suction  surface  at  the  point  of 
tangency  to  the  straight  segment,  (point  I) 
intersects  the  pressure  surface  of  the  adjacent  blade 
at  the  trailing  edge  (point  C),  or  downstream  of 
point  C,  depending  on  the  value  of  the  parameter 

(  C*> 

Referring  to  figure  C2  and  C3  the  following 
relationships  are  obtained: 

xp  =  xc  -  S'  cos  0  C(l) 

YF  =  Yc    *     S'  sin  #  C(2) 

where  &'    is  the  leading  edge  and  trailing  edge  thickness. 

The  value  of  S'  has  been  made  dimensionless  by  dividing  by 
RQ.   All  other  quantities  have  been  made  dimensionless  in 
the  procedures  given  in  appendix  B. 

xt  =  xc  -  S/2  cos  0  C(3) 

yt  =  Yc    *     5/2  sin  0  C(4) 

xe  =  - S  C(5) 

ye  =  Ri  C(6) 

xl  =  -S/2  C(7) 

yL  =  Ri  C(8) 


SS 


The  wedge  angle,  (  is  specified.   With  (  known,  the 
length  of  the  straight  section  <  x   )  from  the  exit  can  be 
found.   For  a  symmetric  blade,  the  length  (.  Jl   )  is  the  same 
at  the  inlet  and  the  exit.   Referring  to  figure  C2,  C3  and 
C4 

J2  =  JL '  <i  -  **>  c<9) 

where  X   =  J(    '    when  the  point  I  lies  on  the  radius  CD' 

(figure  C4>  .   The  slope  of  the  line  C'F  < u)  is  given  by, 

\i    -    arc  tan  <  vf  -  y^')  C(10) 

<  xp  -  xc' > 

The  coordinates  of  C  are  found  by  solving  for  the 
previous  blade  and  rotating  the  coordinates.   These  values 
are  placed  in  an  array  until  needed.   A  straight  line  is 
drawn  from  the  current  blade  to  the  previous  blade  as  is 
defined  as 

C'F  =  <  CyF  -  yc12    *     Cxp  -  xc']2  )  *         C<11> 

Referring  to  figure  C3,  the  slope  of  FG  is  given  by, 

A  =  toe  -  i>    -  £  ]  C(12) 

A  relationship  of  X  '  »    C'F  and  the  line  FG  is  given  by, 

X*  =    C'F  cos  CFG 

=  C'F  cos  C  a  -  y>     -     i    +  u]     C<13) 

From  equation  C(13)  the  coordinates  of  I  and  H  are  found, 

xi  =  xf  -  t    cos  ta  -  ^  -  (]  C(14) 
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yi  =  yF  -  2    sin  ca  -  y>    -   a  c<i5> 

xh  =  -  C  S  +  JL      sin  il  C<16) 

yH  =  Ri/R0  ♦  H     cos  €  C<17) 
The  centers  of  the  suction  side  are  given  by  either, 

XD'  =  xH  *  rz  cos  *  C<18) 

yD'  =  yH  +  R2  sin  £  .  C<19> 
or, 

xq'  =  xi  ♦  Rz  sin  C  a  -  J.  -     £3             C(20) 

yD'  =  yi  -  Rz  cos  c  a  -  y-    -   a  c<2i) 

From  equations  C<17),  C(19)  and  C(ll)  Rz  can  be  found 

since 

XI  "  XH  =  Rz  (cos  C  -  sin  a,  ) 

yielding, 

Rz  =    xt  -  xh  C(22) 

cos  £  -    sin  & 

To  calculate  the  position  of  I  for  a  given  throat 

dimension  (a')  at  the  exit,  further  manipulation  is 

required.   Referring  to  figure  C4,  where  "a"  denotes  the 

center  of  the  exit  throat, 

xa  =  xc'  +  a'  sin  a  C<23) 

and 

ya  =  yc'  "  a'  cos  a  C<24) 

From  equations  C<21),  C<13)  and  C<15), 

Rz  =  <  xf  -  H     cos  A  +     s  +   J?  sin  £  ) 
(cos  £  -    sin  ^) 

=  xf  ♦  S  ■*•    Csin  €  -  cos  A  ]  C<25) 

Ccos  C-sin  ^) 
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Referring  to  figure  C4 
C  R2  ♦  a']2  =  [  XD>-  xa  ]2  +  [  ya  -  yD']2         C(26) 
Using  equation  C(20)  and  C(21)  the  above  equation  becomes, 

[  Rz  +  a'32  =  Cxi-R2  sin  ^  -  xa3 2  *  Cya-yi-R2  cos  AJ 2 
referring  to  equations  C(14),  and  C<15), 

C  R2  ♦  a'] 2  =  [xp  -  X  cos  A  *  Rz  sin  A  -xa]2 

♦  [  ya  -  YF  +  A  sin  ^  -  Rz  cos  A  3  2 
which  reduces  to 

2R2a'  ♦  a'2  =  (xp  -  xa)2  ♦  <ya  -  yf )2  ♦  jj  2        C<27) 
♦2  X  cos  A  <xa  "  xp)  +  2  J  sin  A  (va  -  YF> 
♦  2R2  cos  ^  <ya  -  yF>  +  2R2  sin  A  <XF  ~  xa> • 
Using  equation  C(25)  the  following  results  can  be 
obtained 

a'2  =  (xp-  xa)2  +  (ya-  yF>2  ♦  j£  2  ♦  2  X  cos  A<xa-xp) 

♦  2  X   sin  A  <ya  -yp> 

♦  2  f  xy    •*■     S  +  ■*  Csin  €  -  cos  A  ]  \   sin  A  <xf-xa) 

»     Ccos  (    -  sin  ^  3       / 

♦  2/  xf  ♦  S  ♦  -xCsin  (  -  cos  A  3  \   cos^<ya-yp) 

v     Ccos  i    -    sin  A3        ' 

xf  ♦  S  +  -*Csin  (  -  cos  A] 


*2(-^T 


cos  {    -  sin  A 3 
Ccos  (  -  sin  A.  3a'2  =  C  cos  £  -  sin  &3  *  C<28) 

<XF  -  Xa)2  ♦   (ya  -  yF)2  *JJ2  \ 

♦  2lcoaA(xa  -  xF  >  ♦  2  isinA(ya-yF)/ 

*  2   xf  +  S  ♦  Alain  i    -  cos  A  3_  sinA(xF-xa) 

*  2      xF  *     5  ♦  J?  Csin  (    -  cos  A  ]   cosA(ya-yF) 

*  2      xF  ♦  S  ♦  A Csin  i    -    cos  A  3   a' 
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and  on  combining  terms, 

2 

Ccos  (    -    sin  £>,  3a'   =  C  cos  (  -  sin  A  ]  C<29) 

+  (2  Ccos  (  -  sin  ^  3  cos  A  (xa  -  xf) 

♦2  Ccos  i    -  sin  A3  sin  A <ya  -  yp)  } X 

♦  2  xp  sin  A(xp  -  xa)  +  2  S  sin  A(xp  -xa) 

♦  2   Ccos  €  -  sin  A  3  sin  A  (xf  -  xa> 

♦  2  xp  cos  A(ya  ~yF>  ♦  2  5  cosA<ya  -  yF> 

♦  2   Ccos  €  -  sin  A  3  cos^(ya  -  yp) 

♦  2  xp  a'  ♦  2  S  a'  +  2   a'  Ccos  C  -  sin  A  3 

Ccos  i    -  sin   3  a'2  =  C  cos  i    -  sinA3  $2  COO) 

♦  {  2  Ccos  i    -sin  A3  cos  A  <xa  -  xp) 

♦  2  Ccos  (  -  sin  A3  sin  A <ya-  yF> 

♦  2  Ccos  i    -    sin  A3  sin  A<xp-  xa) 

♦  2  Ccos  £  -  sin  A3  cos  A(ya  -yF> 

♦  2  a'  C  cos  i    -    sin  ^3  ) 

♦  2  xf  sin  A (xp  -  xa>  +28  sin  ACxp  _  xa) 

♦  2  xf  cos  A<ya  "  YF*  ♦  2  S  cos  A(ya  -  yp) 

♦  2    xf  a'  «•  2  S  a' 

The  above  equation  is  quadratic  in  the  unknown  )i  ,    and  a 
solution  can  be  written  explicitly.   The  value  of  Jl    can  be 
found  using  equation  CO)  .   The  result  is  as  follows. 

If 

B    =    {       2    Ccos    i    -    sin  A  3    cosA(xa    -    xp)  COD 
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♦  2  Ccos  (    -  sin  A]  sin  A  <ya  -  yp-) 

♦  2  Ccos  (  -  sin  A3  sin  A<xp  -  xa) 
+  2  Ccos  i    -  sin  A]  cos  A  <ya  -  yF> 

♦  2  a  '  C  cos  i    -  sin  A]  }/  j  'Icos  i    -  sin  A  3 

Let  C  =   2  xf  sin  L  <xp  -xa)  +25  sin  ACxf  -  xa>        C<32) 

♦  2  xp  cos  &<ya  ~YF>  +25  cos  A <ya  -  yF> 

♦  2  xf  a'  +  2  S  a')/  J)'  2  Ccos  i    -  sin   ] 
then  our  quadratic  equation  is, 

O  =  J['  Ccos  £  -  sin  _]  <l-£'>  ♦  2B  JJ '  Ccos  £  -sin   3  •  C<33) 

<1  -  (')  ♦  C  JJ'2Ccos  i    -sin   ]  C<34) 

and  is  of  the  form, 

2 
0=<1-C)   *2B<l-(')+C  C<35) 

so  that 


(l-(')  =   -2B/2  ±  V<2B)2    -  4C 

2 


t  VT 


=   -  B    ±  Y<  B   -  C>  C<36) 

The  following  possibilities  occur: 

1.  If  the  term  under  the  square  root  is  negative  then  a 
solution  does  not  exist. 

2.  If  the  term  under  the  square  root  is  positive,  and 
the  value  of  (1  -  (')  is  between  0  and  1  then  a 
solution  exists',  the  solution  is 


.^r2 


<  1  -  €1  )  =  -B  +  V<B    -  C) 
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3.  If  the  value  of  (1  -  *1>  obtained  in  step  (2)  is 
outside  the  range  of  0  to  1  then  the  solution  is, 


-£-2 


<  1  -  (1)  =  -B  -  1<B  -  C) 
4.  If  the  value  of  (1  - il)    obtained  in  step  (3)  is 

outside  the  range  of  0  to  1  then  no  solution  exists. 

Once  ('  has  been  determined  and  all  other  values  have 

been  found,  the  design  is  complete  and  the  blades  can  be 

drawn. 

Referring  to  figure  C2,  to  draw  the  radius  on  the 

leading  edge  the  points 

x  =  xl  +  S/2  cos  9  C<37) 

y  =  YL  "  S/2  sin  9  C<38) 

are  plotted  over  the  interval  of  9  between  0  and  180 

degrees.   The  radius  on  the  trailing  edge  is  given  by, 

x  =  xt  ~  S/2  cos  ©  C<39) 

y  =  Yf    +  S/2  sin  e  CC40) 

plotted  over  the  interval  of  ©  between  0  and  0+180. 

To  draw  successive  blades  on  the  rotor  the  procedure  is 

to  rotate  the  axes  after  each  blade  is  drawn.   The  procedure 

is  described  by  equation  B<55)  and  B(56). 
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C.2   COMPUTER  PROGRAM  "AWEDG" 

C.2.1   Description 

The  wedge  arc  blade  design  program,  "AWEDG" 
calculates  and  draws  wedge  arc  blading  for  an  initially 
radial-outflow  rotor.   The  difference  between  "AWEDG"  and 
"ACIRC"  is  not  only  the  specified  input  parameters;  "AWEDG" 
calculates  the  straight  section  needed  on  the  suction  side 
to  constrict  the  flow  to  a  know  dimension  at  the  trailing 
edge  tip  of  the  adjacent  blade.   The  program,  for  simplifi- 
cation, is  divided  into  seven  parts  and  each  part  will  be 
described. 

The  first  part  (lines  1  to  460)  contains  the 
values  of  all  specified  parameters.   On  running,  it  also 
allows  the  user  to  specify  the  angles  *P      and  0  to  determine 
early  in  the  program  what  the  value  of  a  will  be.   If  the 
user  wishes  to  change  the  value  of  a  he  can  do  so.   Chord 
length,  camber  radius  and  a  are  calculated  in  this  section. 

The  second  part   (lines  470  to  640)  prints  some 
of  the  output  parameters  and  calculates  the  center 
coordinates  of  the  camber  line. 

In  order  to  calculate  the  outlet  passage  size 
the  preceding  passage  must  be  analyzed.   The  program  does 
this  in  part  three  (line  650  to  800)  by  calculating  the 
position  of  the  trailing  edge  tip  of  the  preceding  pressure 
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line.   Point  C  on  figure  C3  is  calculated  by  rotation  of 
axes. 

In  part  four  (lines  810  to  1020)  the  position  of 
points  C,  F,  T,  E,  and  L,  shown  in  figure  C3,  are 
calculated. 

Part  five  (beginning  on  line  1030)  starts  the 
plotting  routine.   The  number  of  the  blades  requested  to  be 
drawn  are  specified.   The  pressure  side  radius  is  calculated 
and  is  drawn.   Calculations  for  the  straight  line  segments 
are  performed  (lines  1340  to  1570).   Based  on  the  determined 
value  of  the  straight  segment  the  coordinates  of  point  I  are 
dervied  (lines  1630  to  2090). 

In  part  six  the  radius  of  the  suction  side  is 
drawn.  Plotting  for  the  leading  and  trailing  edge  radius 
are  also  accomplished  in  part  six. 

Part  seven  (2100  to  2140)  concludes  the 
procedure  by  printing  the  value  of  the  throat  diameter, 
length  of  the  straight  line  segmnent  and  radii  of  pressure 
and  suction  circle. 

C.2.2    Required  Inputs 

The  following  inputs  are  required: 

Outer  radius 

Inner  radius 

Total  number  of  blades 

angle  p 


96 


angle  0  (must  be  consistent  with^-see 

discussion) 

Diameter  of  the  circle  at  the  passage  exit 

Leading  edge  and  Trailing  Edge  thickness 

Leading  edge  and  Trailing  Edge  wedge  angle 

Number  of  blades  requested  to  be  drawn 

C.2.3    Listed  Outputs 

The  following  quantities  are  output: 

Outer  radius 

Inner  radius 

Number  of  blades 

Angle  ^ 

Angle  0 

Angle  a 

Radial  solidity 

Chord  length 

Thickness  to  chord  ratio 

Leading  edge  and  trailing  edge  thickness 

Leading  edge  and  trailing  edge  wedge  angle 

Throat  width  at  exit 

Radius  of  the  pressure  surface 

Radius  of  the  suction  surface 

Length  of  the  straight  section  on  the 

suction  side. 
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C.2.4    Program  Listing 

Listed  here  is  the  progran  "AWEDG"  for  the 
design  of  wedge  shaped  are  blades.   Table  C2  identifies  the 
symbols  used  in  the  program  with  the  defined  parameters. 
Figure  C6  and  C7  illustrates  computer  generated  shapes  for 
an  outer  radius  of  4  inches,  an  inner  radius  of  3.2  inches 
and  60  blades.   The  value  of  ¥     and  0  are  specified 
parameters  to  limit  the  efflux  angle  a  to  between  S  and  9 
degrees. 


10  REM*******flUEDG  1  AUGUST  1984********** 

29  REM***PROGRfiM  TO  CflLCULRTE  &    DRAW  ARC  BLADING*** 

36  REM***ARCS  AND  STRAIGHT  BACKING  FOR  AN*** 

35  REM***INITIALLY  OUTFLOW  CASCADE***** 

40  DEG 

50    DIN    GC  100 3 i HE  19&]»UC  IOOIjVC  100] 

60    SCALE    -1 j 1 j -it  1 

70  DISP  "SET  PLOTTER  700  BY  7(Y)-<C0NT> "  5 

SO  STOP 

100  DISP  "ENTER  THE  OUTER  RADIUS"; 

110  INPUT  R0 

120  DISP  "ENTER  THE  INNER  RADIUS" 5 

130  INPUT  Rl 

140  R4=R1/R0 

150  DISP  "ENTER  THE  NUMBER  OF  BLADES" 5 

160  INPUT  Z0 

170  L=0 

ISO  S5  =  3t"0/Z0 

190  PI =3. 141592654 

200  S2=2*P1*R0/Z0 

210  R9=(R0-R1 )/S2 

220  S1=(R0-R1>/S2 

230  S0  =  S2/'R0 

240  R2=S0*S1 

250  DISP  "ENTER  PSI-<  BTWN  LE  &  TE  0-360/2"! 

260  INPUT  S3 

270  DISP  "ENTER  PHI -PRESSURE  ANGLE/CANBER< " 5 

2S0  INPUT  TO 

290  A0=90+S3-T0 

300  PRINT  "ALHPA  IS  DESIRED  BWTN  8-9  DEGREES.  THE  VALUE' 
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310  PRINT  "OF  PSI&PHI  CHOSEN  WILL  YIELD  RN" 

320  PRINT  "fiLPHA="R0 "DEGREES" 

330  PRINT  "PSI="S3"AND  PHI="TO 

340  DISP  "WISH  TO  CHANGE  PSI&PHI  Y=1.N=0"5 

350  INPUT  01 

360  IF  01=1  THEN  250 

370  DISP  "ENTER  THE  LE  S,  TE  THICKNESS"? 

330  INPUT  Ql 

390  M0  =  Q1 -RO 

400  DISP  "ENTER  LE  &  TE  WEDGE  ANGLE" 5 

410  INPUT  Ml 

420  C 1  =  R @ * S I N S 3 ••' S  I N  <  T 0 / 2 ') 

430  C0  =  C1-"'R0 

449  r 3 = r:  Q  ..■■•  ( 2  * S I N  (.  T  0  •••'  2  >  > 

450  R4=R1/R9 
460  0  3  =  R  3 * R 0 

47" 0  REM**^***********************^******"*-*"*** 

4S0  REM*******PRINT  STATEMENTS*************** 

490  PRINT  "R0="R0"RI="R1"N0.  BLADES*" Z0 

500  PRINT  "PSI="S3 

510  PRINT  "PHI="TO 

520  PRINT  "CHORD* "CI 

530  PRINT  "ALPHA="R0 

540  PRINT  "RADIAL  30LIDITY="R9 

550  PRINT  "LE  %.    TE  THIC :KNESS=  "Ql 

560  PRINT  "LE  %■.    TE  WEDGE  ANGLE*" Ml 

570  DISP  "ENTER  THE  THROAT  LENGTH  A*"? 

5S0  INPUT  A5 

590  REM******END  PRINT  STATEMENTS***** 

600   REM***^  +  ****  +  ****  +  **^  +  **'*  +  ':,;■'*■*****  +  +  *■'"","'■,■■** 

610  REM****CENTERS  CAN  NOW  BE  DETERMINED**** 

620  X3  =  R3 

630  Y3=R4 

640  REM****DRAW  BLADE  SURFACES**** 

650  REM***FIRST  CALCULATE  COORDINATE  OF  C  ON  BLADE  NO 

660  N9=16 

670  N3=7 

6S0  N1=Z0 

690  1=1 

700  FOR  T4=0  TO  TO  STEP  T0-N9 

7 1 0  X  4  =  ';■  <  3  -  (.  R  3  *  C  0  S  <  T  4  )  ) 

720  Y4=Y3+<R3*SIH<T4>) 

730  X9=X4*C0S < < N 1 - 1 >  *S5  >  +  < Y4*S IN<<N1-1)# S5 > ) 

740  Y9  =  Y4*C0S  (.  ( H 1  -1  > *S5>-  <  X4*S  I N ( <  N 1  -  1 )  *S5 )  > 

750  PLOT  X9,Y9 

760  UC I ]=X9 

770  VC I  ]  =  Y9 

7S0  1=1+1 

790  NEXT  T4 

300  PEN 


Z-l 
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305 
816 

329 
330 
340 
350 
3  6  0 
370 
330 
390 
9  0  0 
910 
920 

9  3  0 
940 
950 
960 
970 
930 
990  r 
1  0  0  0 
1010 

1 0  2  0 
1 0  3  0 
1040 
1050 
1060 
1070 
1030 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1130 
1190 
1200 
1218 
122Q 
1230 
1240 
1250 
1260 
1270 
1238 


REM***F  &  C  ON  PRESSURE  LINE*** 

REM***SECOND  CALCULATED  THE  COORD  I  NATES  OF*** 

N0=1 

1  =  1 

FOR    T4  =  0    TO    TO    STEP    T0--N9 

X4  =  X3-(R3*C0S<T4>  > 

Y4=Y3+(R3*SIN(T4:0 


;5>+<Y4*SIN<<N0-l>*S! 
I5)-<X4*SIN<<N0-1)*S! 


X9=X4*CQS< (NO-1 

Y9=Y4*C0S( (NO-1 

GC I  ]  =  X9 

HE  I  ]  =  Y9 

1  =  1  +  1 

NEXT    T4 

M2  =  GCN9+1  ] 

M3=M2-M0*C0ST0 

M4=M2-<MQ/2>*C0ST0 

L2  =  HCN9+1  ] 

L3=L2+M0*SINT8 

L4=L2+< NO/2 >*S INTO 

M5=-M0 

M6=-(M0/2> 

L5=R1--'R0 

L6  =  R1--R0 

REM*****PL0TTING  SECT  10 

DISP  "INPUT  THE  NUMBER 

INPUT  N 

REM EDIT  NEXT  LINE 

N9=16 

FOR  N0=1  TO  N 

REM****FOR  PRESSURE  LINE**** 

1  =  1 

FOR  T4=0  TO  TO  STEP  T0/N9 

X4=X3-(R3*C0S<T4) ) 

Y4=Y3+<R3*SIN(T4>) 

X9  =  X4*C0S  < < NO- 1 >  *S5  >  +  < Y4  *S I N ( 

Y  9  =  Y 4  *  C  0  S ( ( NO-1 >  *  S 5 ) - < X  4  *  S I N ( 

GC I 3=X9 

HE  I  ]  =  Y9 

1  =  1  +  1 

PLOT  X9»Y9 

NEXT  T4 

PEN 

REM***SUCTION  SIDE 

REM***VALUE  OF  L  3 

REM***TO  GIVE  THE 

E1=L3-VCN9+1  1 

E2=M3-UCN9+1 ] 

L0  =  flTN(-El/'E2> 

C3  =  SQR<E11'2  +  E2t2) 


j ..'  .' 


OF  BLADES  TO  DRAW'S 
FOR  PLOT  RESOLUTION' 


N8-1)*S5>> 

N0-l>*35>;' 


LINE***** 
OLVED  FOR  WILL  BE 
DESIRED  THROAT  AT 


GENERATED*** 
E ;  I  T******* 


100 


298 
300 
310 
3  1 5 
3  2  0 
330 
340 
350 
360 
370 
3  S  0 
390 
400 
410 
420 
430 
440 
450 
460 
470 
430 
490 
500 
510 
520 
530 
540 
550 
560 
570 
530 
590 
6  0  0 
6  1 0 
620 
630 
640 
650 
660 
670 
630 
690 
700 
710 
'20 
'30 
'40 
'50 


) +C0SH 1  *  < V6-L3  > -R4 >  > 


Iil=fi0-S3-Ml 

B9=C3*C0S<  R0-S3+L0-M1 ) 

REM***CALCULATION  OF  BLADE  SUCTION  SURFACE*** 

REM***FQR  A  GIVEN  THRORT  SIZE  RT  EXIT*** 

REM***R4  IS  THE  RADIUS  RT  THE  THRORT  R5-"2*** 

R4=R5/(2*R0> 

X6=UCN9+1  ]  +  R4*SINR0 

Y  6  =  V  C  N  9  + 1  3  -  A  4  *  C  0  S  R  0 

Z5=  < S I ND 1  *  < Y6-L3  >+COSD 1  * ( X6-M3 >  >  *  < COSM 1 -S I HD  1  > 

Z6=  ( S I  HD  1  *  (.  M3-X6  >  +COSD 1  *  <  Y6-L3  >  -R4  )  *  <  S  I NM 1  -COSD 1  > 

Q 1  =  < Z5  +  Z6 )  /  (.  C3*  ( COSM  1  -S I  ND  1 

Z7=  < 2*  <  M3  +  M0  >  *  (.  S I  ND  1  *  ( M3-X6 

Z8=  ( (.  M3-X6 )  t2+  ( Y6-L3 >  T2- A4t2  >  *  ( COSM  1  -S  I  ND  1 ) 

Q2=<Z7+Z3 ) / < C3t2* < COSM 1 -S I ND 1 ) ) 

03  =  0.1  t2-Q2 

IF  Q3<0  THEN  1560 

Q  4  =  -  Q  1  +  S  Q  R  (  Q  3  > 

IF  0<Q4<1  THEN  1470 

IF  Q4<0  OR  Q4>1  THEN  1500 

REM***FIRST  SOLUTION  FOR  E=( 1-Q4>*** 

T9=l-Q4 

GOTO  1550 

Q5  =  -Q1-SQR(C!3> 

IF  0<Q5<1  THEN  1530 

IF  Q5<0  OR  Q5>1  THEN  1560 

REM***SECOND  SOLUTION  FOR  E=(l-Q5>*** 

T9=l-Q5 

GOTO  1570 

PRINT  "NO  SOLUTION  EXISTS" 

BS=B9*<1-T9> 

M7=M3-BS*C0SD1 

L7=L3-B3*SIND1 

M3=-<M0+BS*SINM1) 

L3=R1/R0+B3*C0SM1 

REM***RRDIUS  OF  SUCTION  SIDE*** 

R  7  =  C  M  7  -  M  8 ')  ■■'  ( C  0  S  M 1  -  S I N  D 1  > 

REM***CENTERS  OF  SUCTION  RADIUS*** 

X5=M3+R7*C0SM1 

Y5=L3+R7*SINM1 

REM***PLOT  ROUTINE  TO  DRAW  SUCTION  SIDE*** 

REM***PLOT  RADIUS  ON  LE  TIP*** 

FOR  T4=0  TO  130  STEP  130-N3 

X  4  =  M  6  +  ( M  0  /  2  >  *  C  0  S  T  4 

Y4=L6-(M0/2)*SINT4 

X9=X4*CQS < < NO- 1 > *S5  >  +  ( Y4*S I N C ( NO- 1 ) *S5  >  > 

Y9=Y4*C0S  < < NO- 1 ) *S5  > - < X4*S I N ( ( NO-  1 ) *S5 ) > 

PLOT  X9iY9 

NEXT  T4 
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1760 

1770 

1730 

1790 

1 S  0  0 

1 8 1 0 

1320 

1330 

1840 

1350 

1360 

1870 

1330 

1390 

1900 

1 9 1 0 

1920 

1930 

1940 

1950 

1  9  6  0 
1970 
1930 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2030 

2  0  9  0 
2100 
2110 
2120 
2130 
2140 


X4=M3 

Y4  =  L3 

X9=X4*C0S  ( ( HO- 1 )  *S5 ')  +  ( Y4*S 1 N  ( <  NO- 1 )  *S5 

Y9=Y4*C0S  (.  ( HO-  1  )#S5 )  -  <  X4*S  I N  <  <  NO- 1  >  *S5 

PLOT  X9»Y9 

REM***PLOT  SUCTION  RADIUS*** 

I1=-N1 

I2=2*M1+T0 

I3=<I2-I1).-N9 

FOR  T4=I1  TO  12  STEP  13 

X4=X5-R?*C0ST4 

Y4=Y5+R7*SIHT4 

X9=X4*COS<<N0-l)*S5)+<Y4*SIN<<N0-l)i 

Y9=Y4*C0S  < < NO- 1 )*S5 ) -<X4*S I N ( < NO- 1 > ■ 

PLOT  X9>Y9 

NEXT  T4 

X4  =  M3 

Y4  =  L3 

X9=X4*C0S <  '•  NO- 1  >*S5  >  +  ( Y4*S  I N  ( (. NO- 1  >  ■ 

Y  9  =  Y  4  *  C  0  S  '■  <  N  0  - 1  )  *  S  5  >  -  <  X  4  *  S I N  (  ( \  \  O  - 1  >  ■ 

PLOT  X9»Y9 

REM***PLOT  RADIUS  ON  THE  TE  TIP*** 

I1  =  T0 

I2=T0+130 

I3=(I2-I1>/N3 

FOR  T4=I1  TO  12  STEP  13 

X4=M4-<M0/£>*COST4 

V4=L4+<M0/2>*3INT4 

X9=X4*C0S < ( HO- 1 > *S5  >  +  ( Y4*S I H < < HO- 1 ) 

Y9=Y4*C0S< < HO-1 >*S5 


i  'i 


■z>  j  >  > 
S5>) 


-<X4*SIH((H0-1> 


PLOT  X9iY9 
NEXT  T4 
PEN 

NEXT  NO 
PRINT  "A*="A5 
PRINT  "L="B8 
PRINT  "RADIUS  OF 
PRINT  "RADIUS  OF 
END 


PRESSURE  CIRCLE  ="03 
SUCTION  CIRCLE*" <R7*R8> 
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TABLE  CI.   List  of  Computer  and  Text  Symbols 
for  Wedge-Arc  Blade. 


'AWEDG" 


TEXT 


IDENTIFICATION 


AO 
A4 

A5 

B8 
B9 

CO 

CI 
Dl 
El 
E2 

Gtl] 

HCI] 

I 

II 
12 
13 


2»a' 


i' 


efflux  angle 

throat  dimension  (radius  of 

outlet  circle) 

throat  diameter  (diameter  of 

exit  circle) 

length  of  straight  section 

length  of  straight  section  to 
be  tangent  to  outlet  circle 

c/R0  dimensionless, 
chord  length 

chord  length 

constant  angle 

constant  value 

constant  value 

array  for  current  pressure 
side  x  coordinates 

array  for  current  pressure 
side  y  coordinates 

counter 

Interval  variable 

Interval  variable 

Interval  variable 
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TABLE  CI.   List  of  Computer  and  Text  Symbols  for 
Wedge-Arc  Blade  (con't). 


•AWEDGE" 


TEXT 


IDENTIFICATION 


L2 


Yc 


final  value  on  present 
pressure  line,  y  coordinate 


L3 


yF 


final  value  on  present 
pressure  line  plus  blade 
thickness,  y  coordinate 


L4 


yt 


center  point  of  trailing 
edge,  y  coordinate 


L5 


yE 


Leading  edge  initial  suction 
side,  y  coordinate 


L6 


yi 


point  of  tangency  for 
suction  side,  y  coordinate 


L7 


yi 


point  of  tangency  for 
suction  side,  y  coordinate 


LS 


yH 


initial  point  for  suction 
side  curve,  y  coordinate 


MO 


Leading  and  trailing  edge 
thickness 


Ml 


Leading  and  trailing  edge 
wedge  angle 


M2 


xc 


final  value  on  present 
pressure  line,  x  coordinate 


M3 


xf 


final  value  on  present 
pressure  line  plus  blade 
thickness,  x  coordinate 


M4 


*t 


center  point  of  trailing 
edge,  x  coordinate 


M5 


xe 


Leading  edge  initial  suction 
side,  x  coordinate 
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TABLE  CI.   List  of  Computer  and  Text  Symbols  for 
Wedge-Arc  Blades  (con't). 


•AWEDGE" 


TEXT 


IDENTIFICATION 


M6 
M7 

M8 

N 
NO 

N9 

01 


Q4 
05 
RO 


XL 
XI 

XH 
N 


03 

Rc 

04 

Rs 

PI 

n 

01 

- 

02 

- 

03 

_ 

center  point  of  leading  edge 

point  of  tangency  for  suction 
side,  x  coordinate 

initial  point  for  suction 
side  curve,  x  coordinate 

Number  of  blades  to  draw 

Number  of  blades  to  draw  for 
loop 

Number  of  points  for  interval 
of  arc 

Equals  1  to  change  f-   and  0, 
0  for  no 

output  print  of  Rc 

output  print  of  Rs 

TT 

constant 

constant 

constant,  value  under  square 
root  term  for  solving  for 
(1 

same  as  03 

same  as  03 

outer  radius 
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TABLE  CI.   List  of  Computer  and  Text  Symbols  for 
Wedge-Arc  Blades  (con't). 


EDGE" 

TEXT 

Rl 

Ri 

R2 

e 

R3 

Rc^Ro 

R4 

Ri/Ro 

R6 

Rp/Ro 

SO 

s 

SI 

o*r 

S2 

s/R0 

S3 

r 

S5 

y>* 

TO 

0 

T4 

- 

T9 

- 

X2 

XH 

X3 


X4 


X9 


Xd 


IDENTIFICATION 

inner  radius 

radial  chord,  <RQ  -  Rj.>/R0 

camber  line  radius, 

di mens ion less 

inner  radius,  dimensionless 

pressure  radius, dimensionless 

peripheral  spacing,  2nRo/Z 

radial  solidity 

dimensionless  peripheral 

spacing 

blade  curvature  with  respect 
to  center  of  rotor 
spacing  angle,  360/Z 

camber  angle 

loop  parameter  to  draw  arcs 

value  of  1-  (' 

pressure  surface  center,  x 

coordinate 

camber  line  center,  x 

coordinate 

variable  for  x  values  in 
continuous  loop 

variable  for  x  values  in 
rotational  loop 
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TABLE  CI.   List  of  Computer  and  Text  Symbols  for 
Wedge-Arc  Blades  (con't). 


"AWEDGE"      TEXT  IDENTIFICATION 

Y2         ydp  pressure  center  y  coordinate 

Y3         yd  camber  center  y  coordinate 

Y4         -  variable  for  y  values  in 

continuous  loop 

Y9         -  variables  for  y  values  in 

rotational  loop 

ZO         Z  number  of  blades 
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Figure  CI.  Relaxed  Double  Circular  Arc  Passage 
with  Pressure  Line  Removed. 
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Figure  C2.   Wedge-Arc  Shaped  Blading  Geometry 
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f  aca 


Figure  C3.   Enlarged  View  of  Trailing  Edge 
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radius  a' 


pressure 
surface 


Note:  CD' is  not  normal  to  the  pressure  surface 


Figure  C4 .   Geometrical  Relationship  Between  Throat 
Dimension  and  Suction  Surface  Radius 
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RQ  =  4    Ri  =  3.2   No.  of  Blades  =  60 

Y>     =    9.17        0    =    90.68 
chord     =    0.896198294 
alpha    =    8.49 

radial  solidity  =  1.909859317 
L.E.  £■  T.E.  thickness  =  0.02 
L.E.  &  T.E.  wedge  angle  =  5 
A»  =  0.05225 
L  =  0.037306915 

radius  of  pressure  circle  =  0.629980625 
radius  of  suction  circle   =  0.476827358 


Figure  C5 . 


Computer  Generated  Shapes  for  a  Cascade  of  Wedge 
Arc  Shaped  Baldes  with  60  Blades. 
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